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ABSTRACT 

We apply the Delaunay Tessellation Field Estimator (DTFE) to reconstruct and an- 
alyze the matter distribution and cosmic velocity flows in the Local Universe on the 
basis of the PSCz galaxy survey. The prime objective of this study is the production 
of optimal resolution three-dimensional maps of the volume-weighted velocity and 
density fields throughout the nearby Universe, the basis for a detailed study of the 
structure and dynamics of the cosmic web at each level probed by underlying galaxy 
sample. Fully volume-covering three-dimensional maps of the density and (volume- 
weighted) velocity fields in the cosmic vicinity, out to a distance of 150 h~ 1 Mpc, 
are presented. Based on the Voronoi and Delaunay tessellation defined by the spa- 
tial galaxy sample, DTFE involves the estimate of density values on the basis of the 
volume of the related Delaunay tetrahedra and the subsequent use of the Delaunay 
tessellation as natural multidimensional (linear) interpolation grid for the correspond- 
ing density and velocity fields throughout the sample volume. The linearized model of 
the spatial galaxy distributi on and the correspondin g peculiar velocities of the PSCz 
galaxy sample, produced by iBranchini et alj (|1999t ). forms the input sample for the 
DTFE study. The DTFE maps reproduce the high-density supercluster regions in op- 
timal detail, both their internal structure as well as their elongated or flattened shape. 
The corresponding velocity flows trace the bulk and shear flows marking the region ex- 
tending from the Pisces-Perseus supercluster, via the Local superclusters, towards the 
Hydra-Centaurus and the Shapley concentration. The most outstanding and unique 
feature of the DTFE maps is the sharply defined radial outflow regions in and around 
underdense voids, marking the dynamical importance of voids in the Local Universe. 
The maximum expansion rate of voids defines a sharp cutoff in the DTFE velocity di- 
vergence pdf. We found that on the basis of this cutoff DTFE manages to consistently 
reproduce the value of fi m « 0.35 underlying the linearized velocity dataset. 

Key words: (03.13.6) Methods: statistical, numerical - (12.12.1) Cosmology: large- 
scale structure of the Universe 



1 INTRODUCTION 

The measurement and analysis of the peculiar velocities of 
galaxies forms a major probe of the cosmic structure for- 
mation process. The galaxy velocities reflect the large-scale 
matter flows which according to the gravitational instabil- 
ity scenario for cosmic structure formation go along with the 
formation and emergence of structure in the Universe. They 
are the response to the gradually unfolding underlying inho- 
mogeneities in the cosmic matter distribution. As a result, 
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the comparison of the induced cosmic flows on Megaparsec 
scales with the underlying matter distribution represents a 
potentially powerful instrument to further our insight into 
the dynamics of cosmic structure formation process and in 
inferring the values of a variety of key cosmological param- 
eters. 

Over the past two decades a major effort has been di- 
rected towar ds compiling larg e samples of galaxy pecu liar 
velocities fsee lDeke] (|l994T ) and lStrauss fc Willickl (|l995l ) for 
reviews of the subject). By opening up the window on the 
dynamics of the cosmic structure formation process the anal- 
ysis of these catalogues of galaxy peculiar velocities has lead 
to enormous progress in our understanding of the process. 
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This is particularly true for scales larger than > 10/i _1 Mpc, 
scales on which structure formation still resides in the linear 
phase of development. In particular the Mark III catalog, 
with an effect ive depth sa 60/i~ 1 Mpc, stands as a landmark 
achievement dWillick et al.lll997l ). 

Despite the successes progress has remained limited 
given the sizeable random and systematic errors that be- 
set catalogues of peculiar velocities of galaxies and clusters. 
It allowed the mapping of the cosmic flow field in only a 
rather restricted volume of the nearby Universe. Perhaps 
even more significant is that the low quality of the peculiar 
velocity data, in combination with their inhomogeneous and 
sparse spatial distribution, has prevented the development of 
an equally compelling and complete view of cosmic dynam- 
ics at the higher spatial resolution needed to probe and un- 
derstand the dynamics of nonlinear structures. It is on these 
scales of a few Megaparsec that the structure formation pro- 
cess has progressed to the more advanced quasi-linear stage 
and has left an imprint of genuinely recognizable cosmologi- 
cal structures and patterns. The characteristic flattened and 
filamentary features of the cosmic web, as well as the under- 
dense void regions, are particularly outstanding examples. 
As yet it has not been possible to map the flows in and 
around these emerging structures. 

Progress in the study of large scale motions will not only 
depend on the availability of a substantially larger, better 
defined and considerably more accurate sample of galaxy pe- 
culiar velocities. Equally important will be the use of a more 
sophisticated machinery to handle the discrete, sparse and 
usually inhomogeneously sampled galaxy velocity datasets. 
It does remain a challenge to transform these into prop- 
erly defined velocity maps throughout the surveyed volume. 
It brings to the fore the issues of interpolation of discrete 
data to any location within a given volume as well as the 
issue of the choice for a proper filtering prescription. As 
a consequence of the specific choices and assumptions in- 
volved with a particular method not all information inher- 
ent in the dataset is preserved, and the resulting analysis 
does usually not include all available large and small scale 
properties of the peculiar velocity field. As for the first is- 
sue, the data interpolation, the measured peculiar velocities 
need to be interpolated into regions devoid of data. A con- 
ventional approach is that of interpolation of the discrete 
data set to regular grid locations by means of a fixed inter- 
polation ke£neL_A_jwtabl_e_e2cample i that of the TSC ker- 
nel (see iHocknev fc Eastwood! Il988h . It is straightforward 
to show, however, that this will result in a ma ss-weighted 
quantity (Bcrnard eau fc van de Weygaert] Il996t ). This has 
confused comparisons with analytically derived results, as 
most of these concern volume-weighted quantities. The sec- 
ond aspect is that of the filter choice for the interpolation 
procedure. One may apply a fixed, global or local filter kernel. 
For sensibly defined field values it is important that data get 
filtered over an appropriately large volume. Amongst other, 
shot-noise effects should be suppressed. Conventionally, the 
galaxy peculiar velocities are filtered on a uniform fixed scale 
set by the requirement that also the sparsest sampled regions 
are properly averaged over. The resulting velocity map is one 
in which all the velocity information on scales smaller than 
the filter size has been filtered out. The most widely applied 
example is that of Gaussian filters (e.g. iBertschinger et aU 



1990| ; lFisher et all 19951; iBernardeau et al.l 19951 ; iBaker et ail 



19981 ; IBranchini et al.lll999l ; iDekel et al.lll999h . 

In this study we investigate the performance of the 
Delaunay Tessellation Field Estimator - DTFE - towards 
recovering the spatial density field and velocity flow in 
the nearby Universe. DTFE is a self-adaptive density es- 
timation and multidimensional interpolation scheme which 
doe s not make use of any ar t ificial, user-specified filter- 
ing [SchaaE^TOnde^xSi3 (|2000T ). Currently DTFE is 
the most elaborate tessellation-based filter and interpolation 
scheme available to the field of astrophysics and astronomy 
(see sect. 13. 2| ). It belongs to a wider generic class of Natural 
Neighbour interpolations schemes. As a linear multidimen- 
sional field interpolation scheme DTFE may be regarded as 
the linear first-order equivalent of the higher order Natural 
Neighbo ur algorithms (NN-neighbour) f or spatial interpola- 
tion (see lSibsonlll98ll ; lOkabe et alJl200Ch . 

DTFE and NN-neighbour schemes are based upon 
the use of the Voronoi and De launay tessella t ions of a 
given spat ial discrete point s et l|Voronoil Il908l ; iDelaunavl 
1 1934 ) (see lOkabe et all ( 2000) for an excellent and exten- 
sive overview). Voronoi and Delaunay tessellations epito- 
mize a purely locally defined division of space. Their self- 
adaptive nature concerns both spatial resolution and lo- 
cal geometry, to which they adjust through their defi- 
nition as region of influence amidst their natural neigh- 
bours. DTFE uses the high level of sensitivity of the 
Voronoi/Delaunay tessellation to the local point distribution 
to pro duce an estimate of the loca l density at each sample 
point (van de Weygaert 1991, 1994). For telling illustrations 
see lSchaap fc van de Weygaert! ()2006h - Subsequently, it uses 
the Delaunay tessellation as multidimensional spatial inter- 
polation grid. Unique within the context of NN-neighbour 
techniques. DTFE includes the explicit extension towards 
the self-adaptive determination of density fields on the basis 
of the spatial distribution o f the p oint set itself, defined by 
ISchaap fc van de Weygaertl (J2000) . Also important is that 
its linear nature ascertains its applicability to large data 
sets involving more than a million points. 

We apply DTFE to the study of the spatial mass dis- 
tribution and velocity field in the local Universe on the 
basis of a sample of galaxy p ositions and velocities in- 
ferred from the PSCz catalog l|Branchini et al.l 1 19991 ). A 
uniform, densely sampled and fully volume- covering sam- 
ple of galaxy peculiar velocities in the Local Universe 
does not (yet) exist. We take the alternative of using the 
galaxy velocities and positions inferred from the uniform 
all-sky PSCz galaxy redshift sample. The PSCz galaxy sam- 
ple is the last and best defined galaxy redsh ift selection 
from the IRAS catalog l|Saunders et alj l2000h . This con- 
text guarantees the near uniform and well-defined depth 
and sky coverage essential for a self-consistent reconstruc- 
tion of density and velocity fields within its volume of 
the nearby Universe (Branchini et al.|[l 999; Schmold t et al.l 
1 19991 ; IBranchini et al.ll200ll ; iTeodoro et al-feoOSl ). Although 
this involves a reconstruction on the basis of linear perturba- 
tion theory and the velocity field therefore does not contain 
proper nonlinear components the spatial patterns outlined 
by the galaxies' positions do provide a reasonable spatial 
configuration to test DTFE over its ability to recover pat- 
terns in the velocity field. 

This case study will demonstrate the virtues and po- 
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tential of the DTFE method in the analysis of observational 
data, in particular in the capacity of DTFE to reveal intri- 
cate anisotropic patterns and that of rendering the multi- 
scale character of both density and velocity fields. The ob- 
jective of this study is first and foremost the production 
of maps of the volume-weighted velocity and density fields 
throughout the nearby Universe. The self-adaptive nature 
of the DTFE technique yields an optimal resolution of the 
most interesting features, and this will enable the identifica- 
tion of characteristic patterns in the flow field and of their 
principal sources in the local density field. By default - to 
the DTFE method - this includes shear flows, (super)cluster 
infall and void outflow. By illuminating their relationship 
to the nearby and surrounding large scale structures the 
DTFE procedure will help in increasing our understanding 
of the dynamics of the formati on of the various structural 
components of the cos mic web l|Bond et al.lll996h (also see 
Ivan de Weygaerdl2002h . 

An additional aspect of the presented study is the 
DTFE performance in a non-uniformly sampled galaxy dis- 
tribution. Because the PSCz galaxy catalog on which our 
analysis is based involves a flux-limited selection of galax- 
ies the sampling density gradually decreases outward. While 
DTFE is capable of correcting the density values accordingly 
it also introduces a differential resolution scale along the re- 
constructed density and velocity field maps. While DTFE 
guarantees the preservation of an optimal share of informa- 
tion on the velocity field, this will automatically diminish 
with the outwardly decreasing selection function and the 
corresponding decrease in sampling density. The spectral 
coverage of the reconstructed DTFE fields will therefore be 
seriously affected towards the outer region of the PSCz vol- 
ume, the velocity field reconstructions less so as they are 
dominated by lower frequency components. 

We introduce the PSCz sample in section 2, along with 
a short description of the linearization process to produce 
the sample of galaxy positions and (peculiar) velocities. This 
is followed by an outline of the DTFE procedure for density 
and velocity fields in section 3. A cosmographic description 
of the spatial structures within the PSCz volume of the local 
Universe is contained in section 4, after which in section 5 we 
turn to the velocity flow field within the same volume. The 
issue of the differential spatial resolution of the DTFE den- 
sity and velocity maps is discussed at some length in section 
6, while we describe some four specific features in the Local 
Universe in section 7. Finally, we discuss the spatial distri- 
bution of the velocity divergence field and its relation with 
the density distribution. Given the rather scant information 
on the velocity shear field we briefly treat its appearance in 
section 9. Section 10 concludes our study with a summary 
and a description of prospects and relation to other work. 



2 PSCz: : 

GALAXY POSITIONS AND VELOCITIES 

For the study of the structure and kinematics of the local 
cosmic neighbourhood we base outselves upon the PSCz cat- 
alog. 



2.1 The PSCz sample 

The IRAS-PSCz catalog dSaunders et alj l200(t) is an ex- 
tension of the 1.2-Jy catalog I Fisher et al.lll995T ). It con- 
tains ~ 15 500 galaxies with a flux at 60/^m larger than 
0.6-Jy. For a full description of the catalog, selection criteria 
and the procedures used to avoid ste llar contaminat i on an d 
galactic cirrus, we refer the reader to lSaunders et all (2000). 
For our purposes the most important characteristics of the 
catalog are the large area sampled (~ 84% of the sky), its 
depth with a median redshift of 8 500 km s _1 , and the dense 
sampling (the mean galaxy separation at 10 000 km s" 1 is 
(/} = 1 000 km s _1 ). It implies that PSCz contains most 
of the gravitationally relevant mass concentrations in our 
local cosmic neighbourhood, certainly sufficient to explain 
and study the cosmic flows within a local volume of radius 
~ 120 /i _1 Mpc. 

Because of the flux-limited nature of the PSCz catalog, 
there is a decrease in the objects' sampling as a function of 
distance from the observer. This is quantified by the radial 
selection function of the catalog, ip(r), where the selection 
function is defined as the fraction of the galaxy number den- 
sity that is observed above the flux limit at some distance 
r. To ascertain the proper number density of objects each 
galaxy is weighed by the inverse of the selection function 
(see eqn. For the s election function we us ed the expres- 
sion tp(r) described in lBranchini et all l|l999t ). 

In addition to this selection function correction, further 
necessary corrections are that for the 16% of the missing sky 
devoid of data, due to the cirrus emission and unobserved ar- 
eas, and that for redsh ift distortions. Here we have followed 
iBranchini et all I (|l999f ) in using their spatially reconstructed 
PSCz catalog. 



2.2 PSCz: : from redshift space to physical space 

To translate the redshift space distribution of galaxies in 
the PSCz catalog into galaxy positions and velocities in real 
space we use base our study on a galaxy sample which has 
been processed by t he linear theory for gravitational insta- 
bility (|Peebledll98Ch . 

Our sample has been obtain e d alon g the lines of method 
I described in IBranchini et alj (|l999T ). The method is a 
specific realization of an iterative techniq ue to minimize 
redshift-space distortions (lYahil et al.lll99l|). For a detaile d 
exposition of the method we refer to Branchini et alj |l999). 
Starting with the original redshift distribution as input and 
computing the corresponding peculiar velocity fields, the 
spatial position of each galaxy is corrected accordingly. This 
is repeated in each iterative step until the spatial distribu- 
tion of the galaxies and the inferred peculiar velocities are 
mutually consistent. In order to minimize the uncertainties 
derived from the lack of information on scales larger than the 
PSCz: catalog, the velocity predictions in the PSCz region 
were made in the Local Group frame [LG]. 

In order to translate the distribution of galaxies into 
a gravitational field, and implied velocity field a value for 
the f3 param eter, the ratio between the linear velocity field 
factor /(O) {Peebles 1980]) and the (supposedly) linear bias 
factor b between galaxy and mass density, needs to be as- 
sumed. A value of /3 = 0.5 has been adopted. The relia- 
bility of the corresponding modelled density and velocity 
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has been confirmed in several studies, including density- 
density a nd velocity-velocity comparisons with ot her sur- 
veys (e.g. iBranchini et al.l[200l] : IZaroubi et al.ll2002h . While 
the reconstruction procedure assumes linearity, one may ex- 
pect its validity to ex tend into the early nonlinear stages 
l|Branchini et al.ll 19991 ). At hindsight, this claim may indeed 
be confirmed given the recovery by our DTFE technique 
of well-behaved mildly nonlinear density and velocity diver- 
gence distribution functions (see fig. [6] and sect. I9.4j l. 

An important condition for a consistent linear recon- 
struction of the velocity and density field is that the den- 
sity fluctuations should have a small amplitude. In order to 
ascertain this the gravity field has been smoothed with a 
tophat filter whose radius is 500 km s _1 within the inner 
50 /i _1 Mpc and subsequently increases following the aver- 
age intergalaxy distance (see fig. [2}. An additional implicit 
requirement for the linearization process is that the vorticity 
modes in the peculiar velocity field are to be minimized. For 
the well-defined and representative sample of PSCz catalog 
this is accomplished for R « 5 h~ Mpc. For comparison, 
it may be worthwhile to rea lize that a Gaussi an kernel of 
12 ft -1 Mpc was employed bv lDekel et alj (|1999T ) to obtain a 
succesfull and representative lineariz ation of the measur ed 
velocity field of the Mark III catalog (|Willick et al.lll997l ). 

A complication for the linearization procedure is that 
of the velocity field around high density regions. In the im- 
mediate surroundings of clusters of galaxies, the inflow of 
matter becomes so fast that t hey t urn into triple- valued 
redshift regions. I Romano- Dfazl (|2004l ) addressed this prob- 
lem by considering two different samples. In one sample 
the triple-valued regions were collapsed, in the other they 
were not. Results showed that the implicit smoothing pro- 
cedure of the linearization procedure minimizes this effect. 
Differences between collapsed and uncollapsed samples are 
less than 10% for the bulk flow and velocity shear compo- 
nents, and consistent at the lcr level. Here the velocity model 
which leaves the triple- valued regions uncollapsed has been 
followed. 

The final product of the linearization velocity recon- 
struction procedure is a galaxy catalog containing for each 
individual galaxy its real spatial position and its peculiar ve- 
locity. Our DTFE analysis of the PSCz density and velocity 
field is based upon this density and linearized velocity field. 



3.1 Voronoi and Delaunay tessellations 

The DTFE method is based upon the Voronoi and Delau- 
nay tessellations of the point sample fsee lOkabe et al. I l2000l . 
, and references therein) . A Delaunay tetrahedron is defined 
by the set of four points whose circumscribing sphere does 
not contain any of the other points in the generating set 
l|Delau nav 1934) (triangles in 2D). The Delaunay tessella- 
tion is intimately related to the Voronoi tessellation of the 
point set, they are each others dual. The Voronoi tessellation 
of a point set is the division of space into mutually disjunct 
polyhedra, each polyhedron consisting of the part of space 
closer to the defining point than any of the other points 
l|Voronoil [l90Sl ; lOkabe et all 120001 ). These spatial volume- 
covering divisions of space into mutually disjunct triangular 
(2-D) or tetrahedral (3-D) cells adapt to the local density 
and geometry of the point distribution. 

Astronomical exploitation of Voronoi and Delaunay 
tessellations have as yet remained rather limited, despite 
the fact that a great many studies are based on observa- 
tions which by their nature represent a discretely and ir- 
regularly sampled probe of an underlying smooth distribu- 
tion. Nonetheless, in recent years the self-adaptive virtues 
of Voronoi tessellations have been recognized in a grad- 
ually increasing stream of astronomical studies. Most of 
these seek to use the density sensitivity of Voronoi cell 
volumes and even Delaunay cell volumes to detect den- 
sity peaks like clusters of galaxies amidst a general back- 
ground of galaxies. In astronomy, the first application of 
Voronoi tesse llations as density /inten s ity es timators was 
forwarded by lEbeling fc Wiedenmannl (| 19931 ) , who sought 
to identify X-ray clus ters as overdensities i n X-ray photon 
count s . More recently, Ramella, et alj (120011 ) .iMarinoni et al.l 
l |2002j ). lKim et all l|2002i ) and lLopes et all l|2004l) used sim- 
ilar algorithms to isolate clusters within catalogues of 
galaxy positions produced by galaxy sky or redshift surveys. 
The explicit and additional ability of the more elaborate 
DTFE forma l ism to trace sharp density contrasts impelled 
iBradac et all |2004) to compute the surface density map for 
a galaxy from the projection of the DTFE volume density 
field. The obtained surface density map was used to compute 
the g ravitational len sing pattern around the object, upon 
which ILi et alj (|2006l ) evaluated the method in its ability to 
trace higher-order singularities. 



3 THE DTFE METHOD 

The Delaunay Tessellation Field Estimator 

(DTFE) interpolation method was i ntroduce d by 

Sch aap fc van de Weygaert! (2000) (see also ISchaadlioOrj ). 
for rendering fully volume-covering and volume-weighted 
fields from a discrete set of sampled field values, including 
the density field as defined by the point sampling itself. 
DTFE involv es an extension of the interpolati o n pro cedure 
described by iBernardeau fc van de Weveaertl (l996), who 
used Delaunay tessellations for the specific purpose of esti- 
mating the cosmic velocity divergence field, and showed the 
method's superior performance with respect to conventional 
interpolation procedures. 



3.2 Local spatial adaptivity 

The point in case for its pattern tracing characteristics 
is provid ed by the right-hand panel of F igure [3] (also see 
fig. 4 in ISchaap fc van de Weygaert] |2006| ). showing the 2- 
dimensional Delaunay triangulation for a section along the 
z-supergalactic plane through the galaxy set which we an- 
alyze in this work, the PSCz catalog. The DTFE method 
exploits these virtues and adapts automatically and in an 
entirely natural fashion to changes in the density or the ge- 
ometry of the distribution of sampling points. Instead of 
involving user-defined filters which are based upon artificial 
smoothing kernels the resulting main virtue of DTFE is that 
it is intrinsically self-adaptive. In essence, it involves filter- 
ing kernels which are defined by the local density and geom- 
etry of the point process or object distribution. The Voronoi 
tessellation is used to obtain optimal local estimates of the 
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Figure 1. DTFE velocity interpolation procedure. The top-right panel presents the Delaunay triangulation of the discrete particle 
positions of the velocity field presented at the top-left panel. The central panels show the velocity gradient computed for each triangle. 
The colour scale corresponds to the amplitude of the determinant of the velocity gradient tensor. The right-hand panel depicts the 
3D representation of the gradient surface. The height of each point corresponds to its velocity amplitude. The DTFE velocity field is 
estimated at the grid points indicated by the coloured grid of the bottom-left panel by assuming a linear variation of the velocity field. 
The bottom-right panel presents the outcome of the DTFE velocity interpolation procedure. 



spatial density (see sect. 8.5 lOkabe et al.l l2000). while the 
tetrahedra of its dual Delaunay tessellation are used as mul- 
tidimensional intervals for linear interpolation of the field 
values sampled or estima ted at the location of the sample 
points (|Okabe et alj|200(l ch. 6). 

On the basis of its interpolation characteristics DTFE 
is the first order version of a wider class of tessellation- 
based multidimensional and entirely local interpolation pro- 
cedur es, commonly known as Natural Nei g hbour Interpola - 
tion j Watson! 1 19921; ISambridge et all Il995l ; ISukumarl 1 19981 ) 
(also see Okabe et al.l l200ol . . ch.~6). In a variety of ap- 
plied sciences the concept of natural neighbour interpolation 
and other advanced tessellation applications have already 
found wide application. Particularly successful and note- 
worthy examples may be found in compu ter visualisation 
and surface rendering , geophysics (see e.g. ISambridge et all 
1995l ; ISambridgell999r i and engineering mechanics (|Sukumaj 



1998 



3.3 DTFE: density and velocity fields 

A crucial a spect of the success of the Delaunay int erpolation 
procedure jBernardeau fe van de Weygaertl 1 19961 ) is that it 
reproduces the volume-weighted velocity field, correcting a 
few fundamental biases in estimates of higher order velocity 
field moments. While, often unintentionally, most conven- 
tional interpolation schemes yield the mass-weighted veloc- 
ity, it are the volume-weighted velocity estimates which fig- 
ure in analytical expressions. This is particularly true for ve- 
locity field perturbation analyses within the context of grav- 
itational instability scenarios based on primordial Gaussian 
perturbations. The validity of the Delaunay (and therefore 
also DTFE) velocity field interpolation scheme is perhaps 
most strongly emphasized by its success in reproducing the 
nonlinear velocity divergence distribution function, specifi- 
cally in its detailed agreement with second order perturba- 
tion theory (|Bernardeau fc van de Weygaertll 19961 ). The de- 
rived distribution function appears to be so accurate that it 
enable s the accurate determination of cosmo logical param- 
eters (Bcrnard eau k, van de Weygaertl I1997T 1. This is even 
more interesting as a mildy non-Gaussian velocity diver- 
gence distribution would enable the breaking of the degen- 
eracy between the cosmic matter density Q. a nd the bias 
b between the matter and galaxy distribution dBernardeaul 
ll994l ; lLokas et al.lll995l) . 

An essential and unique step of the DTFE procedure 
concerns the determination of field values at the sample 
points. For the velocity field this simply involves the mea- 
sured velocities at the sample points. More complex is the 
issue for the density/intensity field. Tessellation-based meth- 
ods for estim ating the densit y have been introduced by 
iBrownl 119651 ) and lOrdl \ 19781 ). On the condition that the 
sample points are sampled proportionally to an underly- 
ing density /intensity field, the Voronoi tessellation is used 



to define optimal - local - density estimates (lOkabe et al.l 
l200d. section 8.5). A minor modification was introduced 
bv lSchaap fc van de Weygaertl f2000): in order to assure the 
mass-conserving nature of the DTFE interpolation proce- 
dure the sample point density estimates relate to the con- 
tiguous Voronoi cell (see fig. [3J . 

Amongst the NN-neighbour schemes, DTFE is unique 
in including the tessellation-based density estimates on the 
basis of the sampling point distribution itself. By restricting 
itself to a linear interpolation scheme it remains feasible to 
apply the method to data sets of millions of points. The 
latter is essential for the viability of the technique within a 
cosmological context. 

3.4 DTFE and Cosmic Structure Formation 

Three major characteristics of the Megaparsec scale Uni- 
verse are the (1) the weblike spatial arrangement of galaxies 
and mass into elongated filaments, sheetlike walls and dense 
compact clusters, (2) the existence of large near-empty void 
regions and (3) the hierarchical nature of the mass distri- 
bution, marked by substructure over a wide range of scales 
and densities. According to the standard paradigm of cosmic 
structure form ation, the theory of gravitational instability 
|Peeblesill980l ) this intricate spatial pattern has emerged as a 
result of the gravitational growth of tiny (Gaussian) density 
perturbations and the accompa nying tiny velocity perturba- 
tions i n the primordial Universe. ISchaap fc van de Weygaertl 
( 2006) demonstrate the ability of DTFE to succesfully repro- 
duce and quantify these key aspects of the nonlinear weblike 
cosmic matter distribution in the Megaparsec universe: 

- The hierarchically structured matter distribution is re- 
solved to the smallest possible resolution scale set by the 
particle number density. 

- DTFE retains the morphology cq. shape of the features 
and patterns in the matter distribution. The characteristic 
anisotropic filamentary and planar features of the Cosmic 
web are fully reproduced in the continuous DTFE density 
field. 

- The near-empty voids in the spatial matter distribution 
are reproduced optimally. Both their flat internal density 
distribution as well as their sharp outline and boundary are 
recovered in detail through the interpolation characteristics 
of the DTFE algorithm as well as by its tendency to s uppress 
the shot noise in these sparsely sampled regions (|Schaap| 
I2006D . 

The relation between the cosmic density and cosmic 
velocity field is an important piece of information on the 
dynamics of the cosmic structure formation process. When 
DTFE is applied to N-body simulations of structure forma- 
tion, it does manage to succesfully reproduce the physical 
and spatial correlation between cosmic density and velocity 
fields for highly nonlinear structures. An early analysis of a 
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GIF iV— bod y simulation d id show a few remarkable exam- 
ples (see e.g. Schaadl2006l ) : the velocity flows in and around 
the cores of high-density regions - for as long as it concerns 
single-stream laminar flows - are traced in detail by the 
DTFE density and velocity maps. Within the same simula- 
tions, the voidlike regions were rendered as super-Hubble ex- 
and i ng bubbles, consistent with our view o f void dynamics 



pand ii 

|ickel 1 19841 ; feheth fc van de Weygaert] I2004T). Considerably 
more detailed recent stud i es (ISchaap fc van de Weygaert! 
120031 ; iRomano-Dfazl |2004| ; ISchaad 120061 ) confirmed that 



DTFE manages to trace the correspondence between veloc- 
ity and density field over a large range of scales, resolving 
both small and large scale features of the velocity field. 

Extrapolating this observation, and encouraged by the 
success of Voronoi-based methods in identifying dark ha- 
los in N-body simulations l|Nevrinck et al . 2005) Erad et al.l 
l|2004l ) used DTFE to assess the six-dimensional phase 
space density distribution of dark halos in cosmological N- 
body simulations. While a fully six- dimensional analysis 
may be computationally cumbersome (|Ascasibar fc Binnevl 
120051 ). the splitting of the problem into a separate spatial 
and velocity-space three-dimensional tessellation may in- 
deed hold promise for an innovative analysis of the dynam- 
ical evolution of dark halos. 



do contain some artefacts, of which the triangular imprint 
of its smoothing kernel is the most pronounced one. 



3.6 The DTFE general reconstruction procedure 

Figure [1] describes the various stages of the DTFE field 
reconstruction procedure. It does so with reference to ve- 
locity field reconstruction. For a sp ecificati o n of t he DTFE 
density field procedure we refer t o ISchaad (|2006l ) (also sec 
ISchaap fc van de Weygaertl 120061 ). In summary, it consists 
of the following sequence of steps 

• Point sample 

Defining the spatial distribution of the point sample: 

+ Density field: 
point sample needs to be a general Poisson process of the 
(supposed) underlying density field, i.e. it needs to be an 
unbiased sample of the underlying density field. 

+ General (non-density) field: 
no stringent requirements upon the stochastic representa- 
tiveness of the sampling process will be necessary except 
that the sample volume is adequately covered. In other 
words, the sample points need not form a uniform sample 
of the underlying density field. 



3.5 DTFE: alternatives 

The performance of DTFE has been tested extensively with 
respect to the performance of the rigid grid TSC method and 
the spatially adaptive SPH method. The TSC procedure is 
rigid with respect to the spatial scale as well as the shape 
of the mass distribution. SPH kernels are adaptive with re- 
spect to the local density of points but lack sensitivity to 
the geometry of the mass distribution. 

ISchaap fc van de Weygaert] l|2006l ) present numerous 
tests of the ability of DTFE to resolve spatial substructure 
and to detect anisotropic features. With respect to spatial 
resolution, DTFE manages not only to resolve maps of self- 
similar Soneira-Peebles models with the highest spatial res- 
olution but also to recover the underlying scaling indices. 
Both the TSC and SPH fields fail completely in reproduc- 
ing the proper scaling properties. While TSC does not pro- 
duce any scaling at all, SPH does manage to reproduce scal- 
ing of the density field over a wide range of spatial scales. 
However, it fails fully in recovering the proper scaling in- 
dices. In addition, DTFE is also the only procedure whose 
density field has an autocorrelation function which agrees 
completely, down to the smallest scales, with the two-point 
correlation function of the Soneira-Peebles point distribu- 
tion. 

Also with respect to its ability to recover weblike fea- 
tures in the cosmic matter distribution the virtue of DTFE 
with respect to the TSC and SPH methods becomes more 
evident. While the SPH procedure fares considerably bet- 
ter than TSC its spherical smoothing kernel does introduce 
some artefacts and deficiencies when it gets to resolving the 
finest features. Defined by a fixed number of neighbours, 
SPH tends to smear out the smallest structures, particu- 
larly for anisotropically shaped ones, and results in features 
which occupy a significantly larger volume than the corre- 
sponding galaxy distribution. However, also DTFE patterns 



• Boundary Conditions 

An important issue, wrt the subsequent Delaunay tessel- 
lation computation and the self-consistency of the DTFE 
density and velocity field reconstructions, is that of the 
assumed boundary conditions. These will determine the 
Delaunay and Voronoi cells that overlap the boundary of 
the sample volume. Dependent upon the sample at hand, a 
variety of options exists: 

+ Vacuum boundary conditions: 
outside the sample volume there are no points. This will 
lead to infinitely extended (contiguous) Voronoi cells 
surrouding sample points near the boundary. Evidently, 
these cells cannot be used for DTFE density field esti- 
mates and field interpolations: the volume of the DTFE 
reconstruction is smaller than that of the sample volume. 

+ Periodic boundary conditions: 
the point sample is supposed to be repeated periodically 
in boundary boxes, defining a toroidal topology for the 
sample volume. The resulting Delaunay and Voronoi 
tessellations are also periodic, their total volume exactly 
equal to the sampel volume. While specific periodi c 
tessellation algorithms do exist (jvan de WeygaertHl994T l. 
this is not yet true for most available routines in standard 
libraries. For the analysis of N-body simulations this is 
the most straightforward and preferrable choice. 

+ Buffer conditions: 
the sample volume box is surrounded by a bufferzone 
filled with a synthetic point sample. The density of the 
synthetic buffer point sample should be related to the 
density in the nearby sample volume. The depth of the 
bufferzone depends on the density of the synthetic point 
sample, it should be sufficiently wide for any Delaunay 
or Voronoi cell related to a sample point not to exceed 
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the bufferzone. A clear conditio n for a sufficiently de ep 
bufferzone has been specified by iNevrinck et all (|2005h . 

When involving velocity field analysis, the velocities 
of the buffer points should also follow the velocity 
distribution of the sample and be in accordance with 
the continuity equation. Relevant examples of possible 
choices are: 

- internal: the analyzed sample is a subsample em- 
bedded within a large sample volume, a sufficient number 
of these sample points outside the analysis volume is taken 
along in the DTFE reconstruction. 

- ran dom cloning t echni que: akin to the technique de- 
scribed bv lYahil etatl (|l99ll ). 

- constrained random field: real izations employing 
the existing correlations in the field (j Bcrtschinger 1987; 
Hoff man fc Rib ak 1991; van de Weygaert fc Bertschingerl 
1996). 



• Delaunay Tessellation 

Construction of the Delaunay tessellation from the point 
sample (see fi g. [3^. While we s t ill us e our own Voronoi- 
Delaunay code Ivan de Weygaert! (|l994T ), at present there is 
a score of efficient library routines available. Particularly 
noteworthy is the CGAL initiative, a large library of compu- 
tational geometry routines 1 



• Field values point sample 

Dependent on whether it concerns the densities at the sam- 
ple points or a measured field value there are two options: 

+ General (non-density) field: 
(Sampled) value of field at sample point. 

+ Density field: The density values at the sampled 
points are determined from the corresponding Voronoi 
tessellations. The estimate of the density at each sample 
point is the normalized inverse of the volume of its con- 
tiguous Voronoi cell Wi of each point i. The contiguous 
Voronoi cell of a point i is the union of all Delaunay 
tetrahedra of which point i forms one of the four vertices 
(see fig. [3] 2nd frame, for an illustration). We recognize 
two applicable situations: 

- uniform sampling process: the point sample is an 
unbiased sample of the underlying density field. Typical 
example is that of iV-body simulation particles. For D- 
dimensional space the density estimate is, 



p( Xi ) = (1 + D) 



(i) 



with uii the weight of sample point i, usually we assume 
the same "mass" for each point. 

- systematic non-uniform sampling process: sampling 
density according to specified selection process. The non- 
uniform sampling process is quantified by an a priori 

1 CGAL is a C++ library of algorithms and data structures for Com- 
putational Geometry, see www . cgal . org 



known selection function tp(x). This situation is typical 
for galaxy surveys, V'( x ) ma y encapsulate differences in 
sampling density if)(a, S) as function of sky position (a, 8), 
as well as the radial redshift selection function ^(r) for 
magnitude- or flux-limited surveys. For D-dimensional 
space the density estimate is , 



p( Xi ) = (1 + D) 



V(x,) V(Wi) 



(2) 



• Field Gradient 

Calculation of the field gradient estimate V/| m in each 
D-dimensional Delaunay simplex m (D = 3: tetrahedron; 
D = 2: triangle) by solving the set of linear equations for 
the field values at the positions of the (D + 1) tetrahedron 
vertices, 



fo fi h h 



To ri r 2 r 3 



(3) 



Evidently, linear interpolation for a field / is only meaning- 
ful when the field does not fluctuate strongly. Particularly 
relevant for velocity field reconstructions is that there should 
be no orbit crossing flows within the volume of the Delau- 
nay cell which would involve multiple velocity values at any 
one location. In other words, DTFE velocity field analysis is 
only significant for laminar flows. 

Note that in the case of the sampled field being the ve- 
locity field v we may not only infer the velocity gradient in 
a Delaunay tetrahedron, but also the directly related quan- 
tities such as the velocity divergence, shear and vorticity. 

• Interpolation. 

The final basic step of the DTFE procedure is the field in- 
terpolation. The processing and postprocessing steps involve 
numerous interpolation calculations, for each of the involved 
locations x. Given a location x, the Delaunay tetrahedron 

m in which it is embedded is determined. On the basis of the 
field gradient V/| m the field value is computed by (linear) 
interpolation, 

/(x) = /(xO + V/| m '(x-xO. ( 4 ) 
In principle, higher-order interpolation procedures are also 

possible. Two relevant procedures are: 

- Spline Interpolation 

- Natural Neighbour Interpolation 
Implementation of natural neighbour interpolations is 
presently in progress. 

• Processing. 

Though basically of the same character for practical 
purposes we make a distinction between straightforward 
processing steps concerning the production of images and 
simple smoothing filtering operations on the one hand, and 
more complex postprocessing on the other hand. The latter 
are treated in the next item. Basic to the processing steps is 
the determination of field values following the interpolation 
procedure(s) outlined above. 

Straightforward "first line" field operations are "Image 
reconstruction" and, subsequently, "Smoothing/Filtering". 
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+ Image reconstruction. 
For a set of image points, usually grid points, determine 
the image value: formally the average field value within 
the corresponding gridcell. In practice a few different 
strategies may be followed, dictated by accuracy require- 
ments. These are: 

- Formal geometric approach: integrate over the field 
values within each gridcell. This implies the calculation of 
the intersection of the relevant Delaunay tetrahedra and 
integration of the (linearly) running field values within the 
intersectio. Subsequently the integrands of each Delaunay 
intersection are added and averaged over the gridcell vol- 
ume. 

- Monte Carlo approach: approximate the integral by 
taking the average over a number of (interpolated) field 
values probed at randomly distributed locations within 
the gridcell around an image point. Finally average over 
the obtained field values within a gridcell. 

- Singular interpolation approach: a reasonable and 
usually satisfactory alternative to the formal geometric 
or Monte Carlo approach is the shortcut to limit the field 
value calculation to that at the (grid) location of the image 
point. This offers a reasonable approximation for gridcells 
which are smaller or comparable to that of intersecting 
Delaunay cells, on the condition the field gradient within 
the cell(s) is not too large. 

+ Smoothing and Filtering: 

- Linear filtering of the field /: convolution of the 
field / with a filter function W s (x,y), usually user- 
specified, 

/a(x) = j /(x')W s (x',y)dx' (5) 
• Post-processing. 

The real potential of DTFE fields may be found in 
sophisticated applications, tuned towards uncovering char- 
acteristics of the reconstructed fields. An important aspect 
of this involves the analysis of structures in the density 
field. Some notable examples are: 

+ Advanced filtering operations. Potentially interest- 
i ng applications are those based on the use of wavelets 
(|Martfnez et al.ll2005r ). 

+ Cluster, Filament and Wall detection by mean s of 
the Multiscale Morphology Filter (|Aragon et al.ll2007l ). 

+ Void ident ification on the basis of the cosmic wate r- 
shed algorithm (jPlaten. van de Weveaert fc Jones! [20071 ). 

+ Halo detectio n in N-body simulations 
(|Nevrinck et all 120051 ). 

+ The computation of 2- D surface densities for the 
study of gravitational lensing (|Bradac et al.|[20o3 ). 

In addition, DTFE enables the simultaneous and combined 
analysis of density fields and other relevant physical fields. 
As it allows the simultaneous determination of density and 
velocity fields, it can serve as the basis for studies of the 
dynamics of structure formation in the cosmos. Its ability to 
detect substructure as well as reproduce the morphology of 
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Figure 2. The (mean) PSCz intergalaxy distance as a function 
of distance (to the Local Group), an expression of the diminishing 
sampling density in the flux-limited PSCz sample. 

cosmic features and objects implies DTFE to be suited for 
assessing their dynamics without having to invoke artificial 
filters. 

+ DTFE as basis for the study of the full phase- 
space structure of structures and objects. The phase-space 
structure dark haloes in cosm ological stru c ture f ormation 
scenarios has been studied bv lArad et alj (|2004l ). 



4 PSCZ: TOWARDS THE 

DTFE DENSITY AND VELOCITY FIELD 

Given the determination of the positions and velocities of the 
P SCz galaxies in our sam ple by the linearization procedure 
of iBranchini et all (|l999h . the continuous volume-weighted 
DTFE density and velocity field - and the corresponding 
velocity divergence and shear field - throughout the sam- 
ple volume is computed following the steps outlined in sec- 
tion [XS] 

In our study we analyze and assess the DTFE den- 
sity and velocity maps within an inner region of ra- 
dius 120 h~ Mpc. It allows a more than qual it ative com- 
parison with t h e ana lysis of IBranchini et al.l (|l999h and 
ISchmoldt et alj l|l999h . Out to this radius the sample con- 
tains 10651 galaxies. To include a sufficiently large re- 
gion enclosing the main sources gravitationally contribut- 
ing to the velocity field in the inner 120 /i _1 Mpc we 
make use of t he data within a sphe r ical volume of radiu s 
180 /i _1 Mpc |Branchini et all 1 19991 : iRomano-Diazl 120041 ). 
Out to this radius the PSCz sample still has sufficient struc- 
ture resolution, without shotnoise becoming dominant. The 
number of galaxies in this '"PSCz -source" sample out to 
this radius is 13432. 

An important factor to take into account in the inter- 
pretation and analysis of the resulting DTFE field maps is 
the fact that the inter-galaxy seperation is an increasing 
function of distance. For the PSC2: sample it rises from near 
« 2 /i _1 Mpc within a radius of 20 ft _1 Mpc to ^ 14 /i _1 Mpc 
at the outer radius of 120 /i -1 Mpc (see fig.0. As a result the 
Delaunay tessellation polyhedra will gradually grow in vol- 
ume as a function of distance. Because it are the polyhedral 
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Figure 3. Modelled peculiar velocity field at the galaxy positions (left-hand panel). The plot represents a slice of 2.5 h —1 Mpc thickness 
centered along the z— supergalactic plane. Numbers indicate the major visible structures along the section. 1- Local supercluster, 2- 
Great attractor region, 3- Pavo-Indus-Telescopium complex, 4- Shaplc y supercluster, 5- Coma cluster, 6- Camelopardalis cluster, 7- 
Pisces- Perseus supercluster, 8- Cetus wall, 9- Sculptor void. (From[Branchini et al. 1999). The right-hand panel shows the corresponding 
Delaunay tessellation of the galaxy distribution. The shadowed regions illustrate the "contiguous Voronoi cell" concept for two galaxies, 
marked PI & P2. 



cells which set the local resolution of the interpolated fields 
we need to take into account that the pure DTFE recon- 
struction maps have a diminishing resolution towards larger 
distances. This makes postprocessing a necessity step for a 
quantitative analysis and/or for an objective assessment of 
the resulting density and velocity fields. 



4.1 PSCz: the Delaunay map 

The spatial positions and peculiar velocities of the galaxies 
in the processed PSCz sample are depicted in the left-hand 
panel of Figure[3] The figure shows the galaxies within a slice 
of 2.5 ft _1 Mpc thickness centered along the z— supergalactic 
plane. The Local Group is located at the origin. Velocity 
vectors, indicating the velocity component within the su- 
pergalactic plane, are plotted at the galaxy positions. The 
size of the velocity vectors is normalized with respect to the 
maximum velocity amplitude within the slice. In the map 
the main large scale structures have been labeled. For a de- 
scription of these we refer to section I5TT1 

For reference the righthand frame depicts the Delaunay 
triangulation defined by the projected positions of the sam- 
ple galaxies within the 2.5 /t -1 Mpc thick slice around Super- 
galactic Plane. By showing the meticulous spatial adaptivity 
of the Delaunay tessellation, it forms a good illustration of 
the rationale behind using it as an estimate for local density. 
The figure emphasizes this by means of the two gray areas, 
the contiguous Voronoi cells surrounding a sample point. 
The inverse of these form the DTFE method's local density 



estimate (see eqn. [TJ . The figure also provides a good im- 
pression of the corresponding non-uniform spatial resolution 
of the Delaunay interpolation grid. The spatial resolution of 
the DTFE interpolation grid follows the density of the sam- 
ple points. 



4.2 PSCz: the DTFE interpolation 

Once the spatial position r n and velocity v n of each galaxy 
n in our PSCz sample has been determined, DTFE will 
be applied towards reconstructing density and velocity 
field. throughout the sample volume. 



4.2.1 PSCz DTFE density values 

Evidently, the calculation of the Delaunay tessellation need 
to be done only once, for both the density and velocity fields. 
While the velocity values v n at the locations r n are part of 
the input information, the density values are estimated from 
the volume V(W) n of the surrounding contiguous Voronoi 
cell W n around each cell n (see fig. [3J . For the 3-dimensional 
density maps of the m agnitude-limited PSCz survey we h ave 
the density estimates (jSchaap fc van de Weyga ert 2000) 

^ 4 

?(Xn) = • (6) 

For the PSCz radial s election function ^ vacz we used the 
expression described in Branchini et al l (|l999h . 



DTFE analysis of the PSCz local Universe 11 



4-2.2 DTFE density and velocity gradients 

The value of the density and velocity field gradient in each 
Delaunay tetrahedron is directly and uniquely determined 
from the location r = (a;, y, z) of the four points forming the 
Delaunay tetrahedra's vertices, ro, ri, r 2 and r$, and the 
value of the estimated density and sampled velocities at each 
of these locations, (p ,v ), (joi,vi), (p 2 , v 2 ) and (p3,v 3 ), 



Vp|, 
Vvl. 




(7) 



The four vertices of the Delaunay tetrahedron are both nec- 
essary and sufficient for computing the entire 3x3 veloc- 
ity gradient tensor dvi/dxj. Evidently, the same holds for 
the density gradient dp/dxj. We define the matrix A is de- 
fined on the basis of the vertex distances (Ax n , Ay n , Az n ) 
(n=l,2,3), 
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(8) 



Similarly defining Av n = v n — Vo (n — 1, 2, 3) and Ap n = 
Pn — po (n = 1, 2, 3) it is straightforward to compute directly 
and simultaneously the density field gradient Vp| m and the 
velocity field gradient Vv| m = dvi/dxj in Delaunay tetra- 
hedron m via the inversion, 
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linearization of the input velocities involves a kernel size 
of v5 ft _1 Mpc 2 for the inner 50 /i _1 Mpc, for larger dis- 
tances it gradually increases with the intergalaxy distance 
l(r) as l(r)/y/Z /i _1 Mpc (for the PSCz intergalaxy distance 
see fig. [2}. As a result, the resolution of the velocity field is 
set to a lower limit of </E h~ 1 Mpc. This is sufficient to assure 
the viability of the DTFE velocity field reconstructions. 



5 PSCZ: 

THE DTFE DENSITY FIELD 

To appreciate the power of the DTFE method, the three- 
dimensional map of the DTFE density field in Figure [4] 
provides a telling illustration. For consistency between den- 
sity field and the velocity field, the density field has been 
smoothed with a Gaussian kernel of \/5 /i -1 Mpc. Note 
that for r<^50 h~ 1 Mpc the resolution of both density and 
velocity field increase proportional to the intergalaxy dis- 
ta nce l(r): while the vel ocity field linearization procedure 
of iBranchini et ail l|l999l ) explicitly involve this as effective 
smoothing radius the DTFE adaptive grid resolution au- 
tomatically scales accordingly (see sect. 17.1] ). The depicted 
isosurface level corresponds to structures at 3 times the 
smoothed mean density. On these scales most outstanding 
features correspond to supercluster or void regions of the 
cosmic web. 3 

Many familiar features are concencentrated near the Z- 
supergalactic plane of the 3D density field. For identification 
of the various features we show the density field in this plane 
by means of the contour colour map in fig. [5] The amplitude 
of the corresponding density values can be inferred from 
the colour bar below the map. The superimposed arrows 
indicate the corresponding DTFE reconstructed velocities 
within the supergalactic plane (see sect. [6]). 

It is interesting to c ompare the dens i ty m ap in fig. [S] 
with the map in fig. 5 of IBranchini et al.l (|1999l ). This also 
involves a density field, of the same PSCz volume, smoothed 
with a Gaussian filter with a smoothing scale which spatially 
adapts itself to the local (galaxy) density. While it shares the 
superb spatial resolution with the DTFE map, the imprint 
of the underlying (spherical) Gaussian filter is particularly 
visible in the boundary regions of voids and su perclusters. 
The work bv lSchaap fc van de Weygaertl (|2006l ) shows that 
its lack of geometric adaptivtiy causes structural artefacts 
near boundaries of low-density voids and anisotropic fila- 
ments. 



4-2.3 DTFE: velocity field constraints 

DTFE interpolation of velocity field values is only feasible 
in regions devoid of multistream flows. As soon as there are 
multiple flows - notably in high-density cluster concentra- 
tions or in the highest density realms of the filamentary and 
planar caustics in the cosmic web - the method breaks down 
and cannot be apllied. 

In the study presented here this is particularly so in 
high-density clusters. The complication can be circumvented 
by filtering the velocities over a sufficiently large region, 
imposing an additional resolution constraint on the DTFE 
velocity field. Implicitly this has actually already been ac- 
complished in the linearization procedure of the velocity 
fields preceding the DTFE processing (see sect. [2}. The 



5.1 The DTFE density field: Cosmography 

The Local Group is located at the origin of the maps in 
fig. |3] fig. 2] and in fig. [S] it is embedded in the surrounding 
Local Supercluster region. On the basis of their gravitational 
influence on the surroundings, manifest in the corresponding 
velocity field, most of the main large-scale structures can 

2 We have used the fact that Rq = Rth I V§ jSuto fc Fuiital 
1 19901) 

3 Notice that here, and throughout this study, the density field S 
is in effect the galaxy density field <5 9 . It may be a biased reflection 
of the true matter density field S m , here parameterized by a linear 
bias factor b. 
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Figure 4. 3D reconstructed DTFE PSCz density field. The field has been smoothed with a Gaussian kernel of 5 ft _1 Mpc (for the 
effective smoothing scale see text, sect. 4.2.3). The isosurface represents structures at 3 times the smoothed the mean density. Notice 
the two huge density concentrations around the 2 = plane, the Pisces-Perseus and Cetus wall complex to the right, and the Hydra- 
Centaurus and Pavo-Indus-Telescopium to the left side. A well delineated bridge connects both structures, the Local Supercluster. For 
feature identity see the numbers in fig. \3\ 



be easily recognized in fig. [3] For the purpose of guidance 
and identification we have inserted labels in the map of the 
supergalactic plane galaxy positions in fig. [3] indicating the 
various large-scale structures in our neighbourhood. 

High-density regions (reddish regions) as well as low 
density ones (dark zones) can be easily recognized along the 
isodensity map in fig. [4] One of the immediate observations 
is that the high-density regions, particularly the highly re- 
solved one in the inner region of the 3D map, tend to be flat- 
tened or elongated, a direct consequence of the DTFE ability 
to trace and reproduce the natural shape and anisotropy of 
structural features. 

Two major matter concentrations along the z — 
plane dominate the field (fig. [5}. The complex formed by 
the Pavo-Indus-Telescopium supercluster and the Hydra- 
Centaurus supercluster with its extension toward the Shap- 
ley concentration dominates the left side region of the su- 
pergalactic plane (the orientation of the figure is such that 
the Pavo-Indus-Telescopium supercluster, at [SGX,SGY]~ 
[—40,-10] fr -1 Mpc, is visible at the front lefthand side 
in the 3D image). Towards the other direction of the Lo- 
cal Supercluster we find a similar outstanding mass con- 



centration. Centrally located is the Pisces-Perseus super- 
cluster, extending out in the Cetus wall to the south 
and the Camelopardalis cluster towards the north. The 
Pisces-Perseus supercluster, clearly visible at [SGX,SGY]~ 
[45,-20] ft _1 Mpc, and the Cetus wall, at [SGX,SGY]r; 
[20, —40] ft _1 Mpc and connecting with the barely visi- 
ble Sculptor wall at [SGX,SGY]^ [0, -80] ft _1 Mpc, over- 
shadow the other other structures. The massive matter 
concentrations on either side of the Local Supercluster 
are connected by a thin filamentary bridge, passing near 
the origin of the map and outlining the Local Superclus- 
ter. The filament runs from the Camelopardalis cluster 
([SGX,SGY]« [45,20] h^Upc) towards the Shapley super- 
cluster ([SGX,SGY]^i [-120,70] /i _1 Mpc) and connects the 
Local Supercluster with the Hydra- Centaurus Supercluster. 

The Hydra-Centaurus supercluster, connecting with the 
Local Supercluster on the northwest side, may be largely 
identified with the Great Attractor region (GA). Also the 
Shapley concentration, further out along the northwest axis, 
may be a major contributor to the velocity flows in the Local 
Supercluster. This can be inferred from Fig. [3] which shows 
that the modelled linearized velocity field does not show 
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Figure 5. DTFE density and velocity fields projected along the z— supergalactic plane in a thin slice. The color bar indicates the plotted 
density scale. The density and velocity field have an effective Gaussian smoothing radius of Rq ~ max y/5,l(r)y/E k _1 Mpc, with l(r) 
the intergalaxy separation distance. Density values are scaled according to the bar at the bottom of the frame, while the velocity vectors 
are scaled such that the (black) velocity vector at the bottom left corresponds to 650 km/s. 



any e vidence for a backflow i nto the Great Attractor region 
Branchini et al.l 



Plionis & Kolokotronis 1998; 



999; 



Rowan-Robinson ct al. 2000; 



Basilakos & Plionis 1998). It 



might be the result of a mere artefact of the reconstruction 
procedure while it did reveal itself in the map computed 
by me ans of the alternative procedure of Schmoldt et al.l 
l|l999h . R ecent studies of the dipole in the X-ray cluster dis- 
tribution l|Kocevski fc Ebelindl2006l ) do indeed provide am- 
ple influence for a major dynamical influence of the Shapley 
concentration, be it that the more localized 2MASS sample 



appears to suggest that the tug of the GA re mains the over- 
ridin g influence for the Local Supercluster l|Erdogdu et all 
l2006t) . 

The Great Attractor complex and its extension to- 
wards the Shapley configuration, the Local Supercluster 
and the Pisces-Perseus complex do define a dynamical 
axis of a suggestive quadrupolar configuration. The tug of 
war between the Great Attractor and the Pisces-Perseus 
will be expressed in terms of a strong compressional tidal 
force perpendicular to the axis, responsible for the fila- 
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Figure 6. Density and velocity probability distribution functions (pdf) for PSCz. In the top row the pdfs are plotted in lin-lin, in 
the bottom row in lin-log. The lin-lin plots emphasize the level of non-Gaussianity of the density and velocity fields. The lin-log plot 
emphasizes prominent features in the distribution, in particular the sharply defined low-density edge, and facilitates comparison between 
the velocity divergence pdf and the density pdf. Left panels: Probability Distribution Function (PDF) of the DTFE density field (solid 
line). For comparison we have included the density PDF from of Branchini et al. (1999, dashed line). Righthand panels: the PDF of the 
velocity divergence (see sect. [9}. The DTFE div v distribution is indicated by means of diamonds, the Branchini et al. (1999) one as 
triangles (only in the lin-log plot) as triangles. Note the good mirror-image correspondence between the density and div v plot for the 
DTFE fields (both lin-lin and lin-log). 



mentary geometry a nd the strong velocity she ar near the 
Local Superclnster ILTlie. Yahil fc Jonesl [l98^) , (also see 
Ivan de Wevgaert fc Bertschingerl Il996l ). In all, the map of 
fig. [4] forms a wonderful illu stration of the the oretical frame- 
work of the cosmic web bv lBond et alj (|l996l ), in which the 
large scale tidal field is the web shaping agent. 

Perpendicular to the dynamical axis defined by the 
Great Attractor complex, Local Supercluster and the Pisces- 
Perseus complex, the most outstanding feature in our lo- 
cal cosmic neighbourhood is the Coma cluster. It is visible 
as the only major concentration at the top of the map, at 
[SGX,SGY]« [0, 70] 7i _1 Mpc, embedded within the Coma 
wall of which the maps also allows a glimpse. 

Also highly interesting are the voids visible in the DTFE 
density map in the Local Universe. These can be read- 
ily identified from the 3D map of fig. [4] - and even bet- 



ter from its underdense mirror image - as large empty 
(and in general roundish) regions. The Sculptor void, sur- 
rounded by the Pavo-Indus-Telescopium complex and the 
Sculptor and Cetus walls, is one of the most salient empty 
features along the supergalactic plane. Its impact on the sur- 
rounding velocity field can be easily recognized (see fig- EJ) - 
The Fornax void is located just at the bottom of the plot 
([SGX,SGY]^ [10,-110] ft _1 Mpc). Also the small void lo- 
cated between the Coma cluster and the Hydra-Centaurus 
region is clearly delineated. 



5.2 The DTFE density field: PDF 

Its ability to resolve small and large scale structure with- 
out loss of spatial resolution and the ability to meaningfully 
interpolate the density field in poorly sampled void regions 
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Figure 7. Probability Distribution Functions for the 3 Cartesian velocity components and velocity amplitude. The red PDFs correspond 
to those computed by means of the DTFE technique. The blue PDFs are those from Branchini et al. (1999) 



allows DTFE to recover the density probability distribution 
function, including both the low density and high density 
end of the pdf. Fig. [6] shows the resulting distribution func- 
tion: the solid lines in the lefthand frames (top: lin-lin, bot- 
tom: lin-log) show the pdf of the DTFE density field. The 
intimate connection with the pdf of the corresponding ve- 
locity divergence distribution will be discussed in sect. 19.41 
By assessing the pdf of the density and flow field restricted 
to the inner 50 /i _1 Mpc of the PSCz volume and observing 
it to be largely in agreement with that for the full sample 
seen here we reassured ourselves that effects resulting from 
the heterogeneous spatial resolution (see sec. IT. T p do not 
substantially affect our results. 

On the positive side of the PDF DTFE manages to trace 
the distribution function down to pdf values of 10 -3 , on the 
negative side down to even lower values. One direct observa- 
tion is that the DTFE pdf does not go down to density values 
of -1. The lin-lin plot in the lefthand frame nicely shows the 
mildly non-Gaussian character of the density field, entirely 
according to expectation. The lin-log plot emphasizes the 
meticulous rendering of the low-density regions: on the low 
density side there is a sharp cutoff at a density threshold of 
5 = —0.8, interestingly close to the density of a spherically 



symmetric shell-crossing void in an EdS universe 4 . It forms a 
powerful confirmation of the power of DTFE to reconstruct 
the density distribution in void regions. 

It is particularly informative to compare the DTFE den- 
sity pdf with that of the lin ear density field reconstruction 
from lBranchini et alj (Il999[). who used a rigid grid-bas ed in- 
terpolation scheme (CIC, Tllocknev fc Eastwoodlll98ct l. The 
dashed lines in the lin-lin and lin-log diagram in fig. [6] are 
the corresponding pdf. Several telli ng differences b e tween 
the DTFE and PSCz density map of iBranchini et all (I1999T ) 
can be identified. The gridbased reconstruction used by the 
latter yields significantly lower density values in void regions, 
in the order of 5 m — 1. It is an expression of the inability of 
the rigid gridbased interpolation to recover meaningful den- 
sity values as a result of the sparsity of sample points: when 
no points are found within the grid kernel, the method yields 
a zero density value. On the side of the high density values 
we see that DTFE recovers systematically higher density 
values. This is a consequence of the DTFE ability to trace 



4 The theoretical expectation for a mature and shell-crossing 
void, with a characteristic inverse tophat density profile, is an 



underdensity of A f» -1. + (l./1.7) d 



-0.8. 
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the density field into the most compact and/or highest den- 
sity regions, without smearing these region out as in the case 
of rigid kernel procedures. One may of course argue that it 
would be more appropriate to compare the DTFE results 
with a correspondingly smoothed CIC field. However, here 
we seek to highlight the fact that DTFE does not need any 
such user-specified tuning and does achieve the required fil- 
tering automatically. 



6 PSCZ: 

THE DTFE COSMIC FLOW FIELD 

The DTFE reconstruction of the continuous volume- 
weighted DTFE velocity field of the PSCz catalog pro- 
ceeded along the steps outlined in section 13.61 Because the 
input velocities have been linearized prior to the DTFE pro- 
cessing there is an implicit limit on the resolution of the 
reconstructed flow field rendering additional smoothing un- 
necessary. It also means that one cannot recover features 
in the flow field on a scale smaller than the kernel size 
(y/E /i _1 Mpc for r<50 /i _1 Mpc, oc l(r)/sqrt5 h~ 1 Mpc for 
r > 50 7i _1 Mpc). 

6.0.1 the Supergalactic Plane 

Figure [5] presents the resulting velocity field by means of 
the projected velocity vectors within the Z- supergalactic 
plane, superposed upon the corresponding DTFE density 
contourmaps. The length of the velocity arrows can be in- 
ferred from the arrow in the lower lefthand corner, which 
corresponds to a velocity of 650 km/s. 

The processed DTFE velocity field reveals intricate de- 
tails along the whole volume. The first impression is that 
of the meticulously detailed DTFE flow field, marked by 
sharply and clearly defined flow regions over the whole su- 
pergalactic plane. Large scale bulk flows, distorted flow pat- 
terns such as shear, expansion and contraction modes of the 
velocity field are clear features uncovered by our DTFE tech- 
nique. DTFE recovers clearly outlined patches marked by 
strong bulk flows, regions with characteristic shear flow pat- 
terns around anisotropically shaped supercluster complexes, 
radial inflow towards a few massive clusters and, perhaps 
most outstanding, strong radial outflows from the under- 
dense void regions. 

The map of fig. [5] shows the success of DTFE in con- 
verting a sample of discretely sampled velocities into a 
sensible volume-covering flow field. In particular its abil- 
ity to interpolate over the low-density and thus sparsely 
sampled regions is striking: the voids show up as re- 
gions marked by a near-spherical outflo. By contrast, 
more conventional schemes, such as TSC or SPH (see 
ISchaap fc van de Weygaertl |2006| ;|R omano-Diazll2004r ) . meet 
substantial problems in defining a sensible field reconstruc- 
tion in low density regions without excessive smoothing and 
thus loss of resolution. At the same time, the local nature 
of the DTFE interpolation guarantees a highly resolved flow 
field in high density regions. 

Overall, there is a tight correspondence with the large 
scale structures in the underlying density distribution. 
While the density field shows features down to a scale of 
V5 /i _1 Mpc (within the inner 50 /t -1 Mpc), the patterns in 



the flow field clearly have a substantially larger coherence 
scale, nearly all in excess of 10 h~ Mpc. Of course, a strong 
correlation between density and velocity field is to be ex- 
pected given the artificial origin of the sample velocities as 
they were generated from the galaxy redshift distribution 
through the linearization procedure. Strictly speaking this 
concerns linear features, although it is interesting to see that 
the correspondence remains almost unanimous for mildly 
nonlinear supercluster and void region. The DTFE velocity 
flow sharply follows the elongated ridge of the Pisces- Perseus 
supercluster. In addition we find the DTFE velocity field to 
contain markedly sharp transition regions between void ex- 
pansion and the flows along the body of a supercluster. We 
should bare in mind that this mainly concerns the outline 
of features in the velocity flow, the corresponding velocity 
values remain linear in character. The latter is an artefact 
of the (artificial) linearization origin of the galaxy velocities 
in our sample. 

Massive bulk motions are particularly concentrated 
near and around the massive structure extending from the 
Local Supercluster (center map) towards the Great Attrac- 
tor region and the Shapley concentration. The DTFE map 
nicely renders this pronounced bulk flow towards the Hydra- 
Centaurus region and shows that it dominates the general 
motions at our Local Group and Local Supercluster. The 
most massive and coherent bulk flows in the supergalac- 
tic plane appear to be connected to the Sculptor void and 
the connected void regions (towards the lefthand side of 
the figure). They are the manifestation of the combina- 
tion of gravitational attraction by the heavy matter concen- 
tration of the Pavo-Indus-Telescopium complex, the more 
distant "Hydra-Centaurus-Shapley ridge", and the effective 
push by the Sculptor void region. Conspicuous shear flows 
can be recognized along the ridge defined by the Cetus 
wall towards the Pisces- Perseus supercluster ([SGX,SGY]~ 
[20, —40] h~ Mpc. A similar strong shear flow is seen along 
the extension of the Hydra-Centaurus supercluster towards 
the Shapley concentration. 

The influence of the Coma cluster is beautifully outlined 
by the strong and near perfect radial infall of the surround- 
ing matter, visible at the top-centre of figure [S] Also the 
velocity field near the Perseus cluster, in the Pisces-Perseus 
supercluster region, does contain a strong radial inflow com- 
ponent. 

Perhaps most outstanding are the radial outflow pat- 
terns in and around voids. The intrinsic suppression of shot 
noise effects through the adaptive spatial interpolation pro- 
cedure of DTFE highlights these important components of 
the Megaparsec flow field and emphasizes the dynamical role 
of voids in organizing the matter distribution of the large 
scale Cosmic Web. 

6.1 Velocity Field: PDF 

We have also compared the PDFs of the Cartesian veloc- 
ity components and the total velocity amplitude. The red 
PDFs in figure [7] represent those computed by means of 
the DTFE technique, while the blue ones are those from 
iBranchini et al.1 (119991 ). 

The DTFE velocity pdfs have a healthy Gaussian ap- 
pearance, as we might have expected for a near linear ve- 
locity field. In this respect it is interesting to notice the 
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stark differences between the D TFE distribution funct ions 
and those of the velocity field of iBranchini et all (|l999l ). In 
particular the velocity component in the ^-direction shows a 
marked deviation. We traced the differences back to the low 
density regions, where rigid gridbased methods have consid- 
erable difficulty in defining reasonable density values. On the 
other hand we also note significant differences in the high 
density regions. Gridbased methods are unable to resolve 
these and tend to average the velocities in these regions over 
the volume of a gridcell. This is reflected in the pdf of the 
velocity amplitude: the DTFE pdf testifies of considerably 
higher velocities in the DTFE velocity field. Note that this 
may yield important implications for the analysis of large 
scale bulk flows within the sample's volume, conventional 
methods may yield unjustifiably biased values. 

6.1.1 the Spatial Flow Field 

It remains a rather challenging task to visualize a three- 
dimensional flow field. We chose to show the spatial struc- 
ture of the DTFE rendered velocity flow within the PSCz 
volume within the three mutually perpendicular central 
slices through the centre of the supergalactic plane (ie. the 
Local Group's location). The combination of the velocity 
flow within the 2.5 /i _1 Mpc thick slices along the super- 
galactic X-Y, the X-Z and the Y-Z supergalactic plane pro- 
vides a reasonable impression of the spatial flow field. To 
understand its connection with the underlying matter dis- 
tribution we also depict the density distribution within the 
same planes. 

To appreciate the structure of the PSCz velocity flow 
and to obtain an impression of the density and velocity field 
contributions on different scales we have included the com- 
plementary set of maps in fig [S] and fig. [5] By means of a 
set of 3 x 2 frames, the latter shows the density and velocity 
field filtered on a scale of Rg = 12 /i -1 Mpc, the former on 
FL G = 3.8 /i _1 Mpc. The density field in the X-Y, X-Z and 
Y-Z central slices is represented by greyscale maps in three 
consecutive rows of the lefthand column, the greyscale val- 
ues in all 3 frames are scaled according to the values in the 
bar atop the figure. The corresponding velocity vector maps 
are shown in the righthand column. The velocity vectors 
in the three planes are normalized such that the reference 
velocity vector depicted atop the three frames corresponds 
to 650km/s. The view on a scale of Rg = 12 /i _1 Mpc is 
that of the large (linear) scale contributions to the flow field 
and identifies the agents of these large scale flows. The con- 
siderably more detailed map of fig [8] has the intention of 
identifying the imprint of mildly nonlinear structural fea- 
tures such as superclusters and voids on the corresponding 
velocity flow. 

The combination of fig [5] and fig. [5] is a potentially rich 
source for the study of the structural dynamics in the local 
Universe. Important for the present study is that it empha- 
sizes three key issues of the DTFE velocity field analysis: 

• the spatial adaptivity of DTFE allows the resolution 
of the velocity flow in and around nonlinear features in the 
density field. 

• the spatial adaptivity of the DTFE kernel to the den- 
sity of sample points implies a dependence of the spatial 
resolution of the (raw) DTFE velocity field on the sampling 



density. As a result high resolution maps do involve a rather 
inheterogeneous spatial resolution over the sample volume. 

• the DTFE velocity field has a more homogeneous reso- 
lution than the corresponding density field. This is a result 
of the larger scale character of cosmic velocity fields. 



7 DTFE MAPS: SPATIAL RESOLUTION 

The cosmic matter distribution and cosmic velocity flows 
include contributions over a wide range of scales. One of 
the virtues of DTFE is its large dynamic range, a re sult 
of its adaptivity (see ISchaap fc van de Weygaertl 120061 ) . In 
principal we are therefore equipped for an assessment of the 
structure of the produced DTFE PSCz density and velocity 
maps in terms of their structure at various scales. However, 
before doing so we need to pay attention to the fact that 
for the flux-limited sample involved the DTFE maps do not 
necessarily involve a homogeneous resolution. 

7.1 Heterogeneous vs. Homogeneous resolution 

The cosmographical presentation of the pure DTFE density 
and velocity map in fig. [S] in sec. 15.11 does not include the 
issue of its rather non-uniform, radially declining, spatial 
resolution. This issue becomes prominent in the comparison 
between the maps of fig. [9] and those of fig. [8] While the 
latter have a uniform but low spatial resolution of Rg = 
12 /i _1 Mpc, the higher resolution maps filtered on a scale 
Rg = 3.8 h~ Mpc do lack such uniform resolution. 

By its nature, DTFE resolves regions of a higher sam- 
pling density on a finer scale, as can be readily inferred from 
the structure of the Delaunay grid in fig. [3] As a result the 
spatial resolution diminishes when the intergalaxy separa- 
tion in the sample increases. The increase of the PSCz in- 
tergalaxy separations as a function of distance (fig. [2| thus 
provies an order of magnitude estimate of the effective res- 
olution scale as function of radial distance from the Local 
Group 5 . 

A resolution scale of Rg ~ 3.8 /i _1 Mpc corresponds to 
a radial depth of ~ 30 h _1 Mpc while R G ~ 12 /i _1 Mpc is 
reached at a distance of 120 /i _1 Mpc. By implication the 
maps in fig.|8]are uniform in the inner 30 /i -1 Mpc, but they 
are marked by a gradually diminishing spatial resolution at 
larger distances. It are in particular the density maps which 
show the gradually fading resolution of the 3.8 /i _1 Mpc 
maps most clearly. By contrast, the Rg = 12 /i _1 Mpc maps 
show a uniformly resolved view of the structure and flow in 
the Local Universe. A scale of Rg = 12 ft -1 Mpc corresponds 
to the sample intergalaxy distance at the outer edge of the 
PSCz sample. 

For a proper quantitative analysis and assessment of the 
DTFE density and velocity field reconstructions it is of the 
utmost importance to take into account the gradually di- 
minishing spatial resolution of the raw DTFE maps. In this 
respect we also need to take into account that cosmic den- 
sity and velocity fields are composed of contributions from 

5 Note that the resolution scale is not exactly a function of radial 
distance, it is entirely determined by the local sampling density. 
On average this does increase as function of radial distance. See 

fig. E] 
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Figure 8. The DTFE density (lefthand) and velocity vector (righthand) maps of the PSCz sample in three mutually perpendicular 
central supergalactic planes: the X-Y, X-Z and Y-Z plane (from top to bottom). The density and velocity field is smoothed on a scale of 
3.8 h _1 Mpc. This implies a uniform resolution up to a radius of 30 h~ 1 Mpc around the centre, with a gradually diminishing resolution 
towards the outer edge of the sample volume. The greyscale bar at the top lefthand indicates the density values. The velocity vectors 
are scaled such that the vector atop the figure corresponds to 650 km/s. 



Figure 9. The DTFE density (lefthand) and velocity vector (righthand) maps in three mutually perpendicular central supergalactic 
planes: the X-Y, X-Z and Y-Z plane (from top to bottom). The density and velocity field is smoothed on a scale of 12 h~ iMpc, involving 
a uniform resolution over the sample volume. The greyscale bar at the top lefthand indicates the density values. The velocity vectors are 
scaled such that the vector atop the figure corresponds to 650 km/s. 



a range of scales. By implicitly filtering out small-scale con- 
tributions in regions with a lower sampling density DTFE 
will include a smaller range of spatial scales contributing to 
a field reconstruction. On average, they will therefore cor- 
respond to regions with a lower amplitude, a subtle point 
of importance in any quantitative analysis. A comparison 
of DTFE density and velocity maps in the PSCz volume 
shows that the latter seems to be hardly affected by this 
effect, a manifestation of the fact that the cosmic velocity 
field is dominated by larger scale modes than the density 
field (|Peebleslll980| ). 



7.1.1 The Large Scale Local Universe: 
the 12 ft _1 Mpc maps 

In the large scale Rg = 12 /i _1 Mpc density and velocity 
field of fig. we recognize nearly all features of the previ- 
ous discussion on the structure of the supergalactic plane 
(sect. 15. ip . The X-Y supergalactic plane clearly involves a 
local matter concentration: while the X-Y plane represents 
a more or less even distribution of in particular high-density 
features, we observe a concentration of high-density regions 
towards the corresponding Z=0 axis in the X-Z density map 
and a more moderate density distribution with some occa- 
sional high-density patches in the Y-Z map. One of the most 
conspicuous mass concentrations in the latter is the Coma 
cluster, at around [SGY,SGZ]^ [70,0] /i _1 Mpc, embedded 
within a massive "Great Wall" complex extending along the 
Z axis. It clearly emphasizes the moderate density values 
near our own LOcal Supercluster, lying in between the mas- 
sive Hydra Centaurus-Shapley concentration and the Pisces- 
Perseus concentration in the X-Y plane. The only massive 
feature which had not been traced in the X-Y map is the 
concentration near [SGY,SGZ]r; [-40,-100] ft _1 Mpc. 

The presence of massive matter concentrations in the 
Rg — 12 h _1 Mpc map induces strongly correlated patterns 
in the flow field. The X-Z map is dominated by a strong 
and rather coherent flow towards the Z= (X-Y) super- 
galactic plane. Interesting is to see that this flow is tied 
in with an almost equally prominent outflow from the low- 
density regions towards the upper righthand of the maps. 
A similar pattern can be recognized in the Y-Z map, where 
near the centre we observe an equally massive flow towards 
the supergalactic Z=0 plane. It appears to be coupled to 
a truely massive radial outflow from the underdensity near 
[SGY,SGZ]^ [0, 80] /i _1 Mpc. If anything, this demonstrates 
the key role of low density regions in the dynamics of the 
local Universe, and in the formation of cosmic structure in 
general. 



Within the X-Y supergalactic plane we recognize 
the familiar flow patterns discussed in the previous sec- 
tion, sect. 15.11 The coherent bulk flow along the Local 
Supercluster-Centaurus-Shapley axis. It is coupled to coher- 
ent towards the Pisces-Perseus region, via a saddle point 
near the Local Group location: the dominant shear flow in 
the Local Universe. Radial outflow low-density regions are 
somewhat less prominent, the outflow from the local Sculp- 
tor void appears to be merged into a large scale coherent 
bulk flow. 



7.1.2 Flow Field scale dependence: 
3.8 /i _1 Mpc vs. 12.0 /i _1 Mpc 

The impact of DTFE is particularly emphasized by the com- 
parison of the medium resolution 3.8 h~ 1 Mpc maps and the 
12.0 /i _1 Mpc maps. The more detailed density map of the 
maps in fig. [8] reveals considerably more detail in both den- 
sity and velocity maps. 

The structure in the X-Y density map has already been 
discussed in detail in sect. 15.11 For the Rg = 3.8 /i -1 Mpc 
density there are two major points of notice in compari- 
son to the Rg = 12 /i _1 Mpc maps. The first is the grad- 
ually decreasing density values towards the outer regions: 
most prominent features are concentrated within the in- 
ner 30 - 40 /i _1 Mpc radius. This is a direct consequence of 
the decreasing spatial resolution of the non-uniform Rg = 
3.8 ft _1 Mpc maps (see sect. |7J. Another feature of inter- 
est is the more prominent presence of (deeper) void regions 
wrt. the corresponding Rg = 12 /i _1 Mpc map. This is in 
particular true for the Y-Z plane. 

The Rg = 3.8 /i _1 Mpc velocity field stands out by con- 
siderably more sharply defined and outlined features. All 
features that were recognized in the Rg = 12 /i -1 Mpc ve- 
locity field can also be identified, be it they are much more 
sharply defined. Moreover, all over the volume these sharp 
velocity features are clearly correlated with features in the 
density field. The mildly nonlinear nature of supercluster 
complexes is reflected in pronounced bulk and shear flows. 
In the velocity map these often go along with remarkably 
sharply defined transition regions where the flow changes 
abruptly as it encounters a flow from another direction. Also 
void outflows appear to be much more prominent than the 
weakly positive divergence regions in the large scale map. 
The DTFE void outflows have the appearance of superhub- 
ble bubbles, and relate very we ll to the deep nonlinear po- 
tenti al wells of true voids (see ISchaap fc van de Weygaert! 
I2007D . In addition, we can recognize small-scale nonlinear 
features that were not seen in the Rg = 12 /i -1 Mpc ve- 
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locity map. The radial infall towards Coma and the radial 
outflow out of the Sculptor void are telling examples. 

Both aspects strongly underline the ability of DTFE 
to trace velocity flows into nonlinear regions. This forms a 
strong argument for DTFE's potential in the analysis of the 
dynamics of mildly nonlinear structures, such as filamentary 
or wall-shaped superclusters an d voids. The DTFE analysis 
of Vir go CDM simulations by ISchaap fc van de Weygaerj 
120071 ) does indeed provide quantitative evidence for this 
ability of DTFE. Only when the nonlinear evolution of a 
structure has proceeded towards a higher nonlinear stage, 
marked by multistream regions and shell crossing displace- 
ments, the linear interpolation scheme of DTFE will break 
down. 



8 INDIVIDUAL COSMIC STRUCTURES 

The ability of the DTFE method to resolve and identify 
both small and large scale features, on the condition of being 
sufficiently sampled, makes it into a highly promising basis 
for the development of structure detection algorithms. 

To illustrate the variety of structures that may be iden- 
tified by DTFE we zoom into four different regions of the 
DTFE PSCz map, each of a different character. The four 
regions are the Local Supercluster, the Coma cluster, the 
Sculptor void and the Pisces-Perseus & Cetus supercluster 
filament. The four corresponding panels in Figure [TU] show 
the density field and corresponding velocity field vectors. 
To ascertain optimal resolution the density field has been 
smoothed with a Gaussian kernel of a mere 1 h~ 1 Mpc. It is 
important to realize that as a result of the irregular spatial 
resolution this does not necessarily mean that all structures 
of that size are revealed. Note that the physical size of the 
corresponding panels varies with the spatial extent of the 
depicted structure. 



8.1 the Local Supercluster 

The top-left panel shows the complex velocity field in and 
around the Local Supercluster. The Virgo cluster is located 
just above the Local Group. Particularly prominent is the 
presence of the Great Attractor region and its dominant in- 
fluence on the velocity field. At the location of the Local 
Group ([SGY,SGZ]« [0,0] ft^Mpc) we find a large bulk 
flow, part of a large scale shear flow towards the Great At- 
tractor region. The shear patterns is a result of the exerted 
gravitational influence of the Great Attractor and Pisces- 
Perseus supercluster citep[e.g. see] [Jemiliophd. 



8.2 the Pisces-Perseus supercluster 

The bottom-righthand frame focusses on the Pisces-Perseus 
supercluster and the Cetus wall. The prominent flow along 
the Pisces-Perseus chain towards the Cetus wall is a clear 
indication of the dynamical connection between these two 
structures. Their gravitational influence can be traced along 
the whole zoomed region. The velocity field around the un- 
derdense region located at [SGX,SGY]« [45, -60] h~ 1 Mpc 
is strongly distorted by these two massive structures. Note 
the presence of the shear pattern near the top-left corner, 
near the location of the Local Group. 



8.3 the Coma Cluster 

The top-right hand panel shows how the Coma cluster, em- 
bedded within the Coma wall, distorts the velocity field in 
its surroundings. The DTFE reconstruction clearly depicts 
the almost isotropic infall into the Coma cluster. The slight 
offset is an artefact of our implementation of the DTFE 
interpolation towards the image grid locations. Coma is em- 
bedded within the Coma wall. The shear pattern visible near 
the bottom-left is the result of the opposing forces between 
the Coma region and the Local Supercluster region. 



8.4 the Sculptor Void 

The Sculptor void and surroundings is the subject of the 
bottom lefthand frame. The lowest measured DTFE den- 
sity contrast value (smoothed at 1 fa -1 Mpc) in this region is 
~ —0.78 at the deepest of the void. At the smoothed scale of 
Vf> /i _1 Mpc the DTFE density t hreshold is -0.74, in agree- 
ment with the reported value bv lPlionis fc Basilakosl {2002) 
of —0.69. The velocity field of this almost "empty" expand- 
ing region is distorted by the surrounding matter distribu- 
tion. Small, yet detectable, distortions delineate the bound- 
ary of this void with the surrounding matter distribution, 
cq. the Sculptor wall. 



9 THE DTFE PSCz VELOCITY DIVERGENCE 

A direct spin-off of the DTFE velocity field analysis is that 
it implicitly comes along with the velocity field divergence, 
shear and vorticity. In the outline of the DTFE method, in 
sect. 13.61 we described how the DTFE method determines 
the gradient V/| m of the sample field / in each Delaunay 
simplex m (eqn.0. For the three spatial components of the 
velocity field, (v x ,v y ,v z ), this translates into a velocity gra- 
dient counting nine components dvi/dxj (eqn. [10}. Within 
the context of DTFE each velocity gradient dvi/dxj has a 
uniform value within each Delaunay tetrahedron. By impli- 
cation, the reconstructed DTFE velocity gradient field is not 
continuous. 



9.1 Velocity divergence, shear and vorticity 

From the 9 velocity gradient components dvi/dxj we can 
directly determine the three velocity deformation modes, the 
velocity divergence V • v, the shear Oij and the vorticity u), 
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where uj — V x v = e 3 Wy (and is the completely anti- 
symmetric tensor). In the theory of gravitational instability, 
there will be no vorticity contribution as long as there has 
not been shell crossing (ie. in the linear and mildly non- 
linear regime). For the evolution of density perturbations, 
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Figure 10. Density and velocity zooms for four different regions along the supergalactic plane indicated by the labels at the top of each 
frame. The density field has been convolved with a Gaussian kernel of 1 fc _1 Mpc for a better impression of such field. Note that the 
density colour contour values in each of the frames are chosen differently, their scaling reflected in the color bars at the right side of each 
frame. The velocity vectors at the top of each frame, at the righthand side of the frame's name, correspond to a velocity of 650 km/s. 



the velocity divergence is an important velocity field com- 
ponent from the early linear regime onward. As nonlinearity 
sets in the deforming character of accompnaying shear flows 
become increasingly decisive in shaping the cosmic matter 
distribution. It i s in particular via the Zel'dovich formalism 
l|Zeldovichlll970h and the Cosmic Web theory by Bond and 
collaborators that we have come to appreciate the key role 
of tides and velocity shear in the formati on of the observe d 
weblike Megaparsec matter distribution l)Bond et al.l fl996). 

The relationship between the cosmic density field on 
the one hand and the velocity divergence and shear fields on 



the other hand contains substantial information on the dy- 
namics and formation of structure in the Universe. Within 
the context of the corresponding DTFE analysis it is im- 
portant to realize that there is an important and implicit 
discrepancy: the DTFE density and velocity fields are en- 
tirely continuous. The velocity divergence fields consist of 
discontinuous patches: the Delaunay tetrahedra, each with 
a constant value of V ■ v. This renders a perfect 1-1 corre- 
spondence between the raw DTFE density and V • v fields 
is unfeasible. Filtering over the scale of a Delaunay simplex 
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Figure 11. Lefthand: DTFE velocity-divergence field, on a smoothing scale of Rq = s/5 h — 1 Mpc projected along the Z-supergalactic 
plane. The thin slice corresponds to the one presented in Fig. [5] The velocity-divergence is in units of the Hubble parameter Hq. The 
color bar represents the scale of the plotted velocity-divergence. Righthand: The map, along the Z-supergalactic plane, of the difference 
between the DTFE velocity-divergence field and the linear theory prediction for the corresponding DTFE density field. The color bar 
indicates the corresponding field values, in units of km s _1 Mpc - 1 . 



Figure 12. The DTFE velocity-divergence field on a smoothing scale of 3.8 h —1 Mpc (lefthand column) and 10 h _1 Mpc (righthand 
column). The maps in the consecutive rows concern the three mutually perpendicular central supergalactic planes: the X-Y, X-Z and Y-Z 
plane (from top to bottom). Note that the map on a scale of 3.8 /i _1 Mpc is uniform only out to a radius of 30 h _1 Mpc, with a gradually 
diminishing resolution towards the outer edge of the sample volume. Also compare the images with the density fields in fig. fig. [8] and [9] 
they are approximately the negatives of the corresponding density fields. The greyscale values are in units of the Hubble parameter Ho, 
with their values given in the bar at the top of the figure. 



is usually necessary to restore the physically expected rela- 
tions. 

In the following subsections we present the maps of the 
velocity divergence field and the velocity shear field in the 
same central supergalactic planes as in figs. [5] and [8] By 
smoothing the velocity divergence, velocity shear and vor- 
ticity fields with the Gaussian kernel Re = v5 ft _1 Mpc, 
akin to the density map in fig. [5] and fig. [9] we restore a 
sense of continuity to the resulting maps. 



9.2 The DTFE velocity divergence map 

The DTFE normalized velocity divergence estimate 6 is the 
sum of the trace of the DTFE velocity-gradient components, 

? s ^x= i f^+^+^y (ID 

H Q Ha \ ax ay oz I 

with Ho the Hubble constant 6 . In Figure QT] (lefthand) we 
have plotted 9. 

The expanding and contracting modes of the velocity 
field are clearly delineated. Expanding regions, correspond- 
ing with or surrounding large and deep underdense voids, 
are identified as those with red to yellow tones marking 
positive divergence modes. The Sculptor and Fornax voids 
and the underdense regions around the Coma cluster can be 
immediately recognized. Other conspicuous expanding re- 
gions are that near the Camelopardalis cluster, located at 
[SGX,SGY]r; [40,45] /i _1 Mpc, and the one at the right of 
the Cetus wall at [SGX,SGY]^ [50,-55] /i _1 Mpc. 

The regions with a negative divergence are contract- 
ing, matter is falling in along one or more directions. The 
strongest contractions, represented by the blue tones, are 
related to the peaks in the density field which can be im- 
mediately identified with the most massive structures lo- 
cated along this slice. The Hydra-Centaurus supercluster is 
the most prominent region of infall. Easily recognizable are 
all clusters such as the Virgo, Camelopardalis and Coma 
cluster, the Pisces-Perseus supercluster, the Cetus wall, the 
Pavo-Indus-Telescopium complex, and partly also the Shap- 
ley concentration. These high-contrast blue regions are em- 
bedded within regions with more moderate infall motions. 
These regions, with green contour values, have densities 



6 We have adopted a value of h = 0.7 (Hq 
100 h km s" 1 Mpc -1 ) 



slightly in excess of the mean. They outline a percolating 
region with a roughly filamentary shape. 

In the interpretation of the velocity divergence map it is 
crucial to take account of an important artefact: high ampli- 
tudes of the velocity divergence are seemingly concentrated 
in the inner region of the PSCz volume. This can be under- 
stood on the same grounds as a similar tendency for the den- 
sity field, the lack of small-scale power in the outer regions 
as a result of the diminishing sampling density. A correction 
for this effect is only possible on the basis of a cosmological 
model assumption and does not fall within the context of 
this study. 



9.3 The density-velocity divergence relation 

The velocity divergence and t he density cont rast are related 
via the continuity equation (Peebles 1980). In the linear 
regime this is a strictly linear one-to-one relation, 

iv.v(x,t) = -^ fl (i){(x,i), (12) 

linking the (galaxy) density perturbation fiel d 8 to the pecu - 
liar velocity field v via the factor /(O m ) (see |Peebleslll980l ) . 
In this we take into account that our discussion concerns 
the density field as sampled by galaxies, whose distribution 
is supposed to be linearly biased with respect to the under- 
lying matter density. There remains a 1-1 relation between 
velocity divergence and density into the mildly nonlinear 
regime (see eqn. I13|l . which explains why the map of the 
velocity divergence is an almost near perfect negative 
image of the density map (see fig. [SJ). 

The approximate validity of the linear divergence- 
density relation (|12[1 . on a scale of Rg = v5 /i~ 1 Mpc, over 
the PSCz volume can be appreciated from the contour map 
in the righthand frame of fig. 1111 The map, in the same su- 
pergalactic X-Y plane as in the lefthand frame, shows the 
difference between the divergence term on the lefthand side 
of eqn. [12] and the density term on the righthand side, in 
units of km s _1 Mpc -1 . Evidently, over the entire volume 
the residual map is near uniform, with the exception of a 
"hotspot" near [SGX,SGY]^ [-10,-10] /i _1 Mpc. If any- 
thing, we see a slight tendency for the differences to increase 
outward. This is most likely to be understood from a differ- 
ent reaction of the velocity field to the diminishing sampling 
density than the density field. 

The residual map's fig. [11] impression of near linearity 
is confirmed by the scatter plot between velocity divergence 
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Figure 13. Scatter plot of density versus velocity divergence, 
in terms of a contour plot of the density of scatter points. The 
solid line indicates the linear density-velocity divergence relation. 
Two approximate nonlinear relations are also indicated, the short- 
dashed line depicts the approximation of Bernardeau (1992), the 
long-dashed one the approximation by Nusser et ai.(1991). 



(V-v) and (galaxy) density 5 in Figure [13] The scatter point 
density in the figure is shown by means of a contour map of 
the scatter point density. 

The scatter plot is compared with three different re- 
lations. Superimposed as a solid line is the linear relation, 
V • v = —Ho/36. The near perfect linearity of the (V • v)-<5 
relation reflects the origin of the velocity estimates on the 
basis of a linearization process. Deviations from linearity, 
be it minor, are observed only for the lowest and highest 
density values. 

Even in the quasi-linear and mildly nonlinear regime 
the one-to-one correspondance between velocity divergence 
and density re mains intact, be it that it involves higher or- 
der terms (see iBernardeau et al.l 12002. for an extensive re- 
view). With i n the context of Eulerian perturbation theory 
IBernardeau! (|l992T l (B) derived an accurate 2nd order ap- 
proximation form the relation between the diverg ence and 
the density perturbation S(x). iNusser et all (|l991f ) (N) de- 
rived a similar quasi-linear approximation within the context 
of the Lagrangian Zel'dovich approximation. According to 
these approximate nonlinear relations, 



'|/(0») [l-(l + «5(x)) 2 / 3 ] (B) 



- V ■ v(x) = { 



-/(a, 



l + 0.18<5(x) 



(13) 
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for a Universe with Hubble parameter H(i) and matter den- 
sity parameter £l m . For comparison we have also included 
these approximate relations for mildly nonlinear density and 
velocity field in fig.[T3](B: the short-dashed line, N: the long- 
dashed line). These relations are clearly too pronounced for 
the velocity divergence map of fig. 1111 



9.4 The Velocity Divergence PDF 

The righthand panels of figure [5] shows the PDF for 
the DTFE velocity divergence field. The nice mildly non- 
Gaussian character of the pdf (in the top lin-lin plot) ar- 



gues for the quality of the DTFE velocity field reconstruc- 
tions: the DTFE velocity field interpolation scheme is able 
to recover to considerable accuracy the velocity divergence 
pdf into the nonlinear regime from a discrete particle dis- 
tribution. It may also be noted that the nice mildly nonlin- 
ear character of the DTFE density and velocity divergence 
field reconstructions are, at hindsight, an i ndication for the 
validit y of the "linearization" procedure of iBranchini et all 
( 1999) into the early nonlinear regime. 

The velocity divergence pdf forms a near perfect mir- 
ror image of the density field pdf (lefthand panel), to be 
expected on the basis of the continuity equation through 
which the velocity divergence and the density field are re- 
lated (see next section |9~3)) . This relation is reflected in the 
density and div v pdfs in the lin-lin plots of the top row and 
is even more clear when comparing the lin-log equivalents in 
the bottom row. DTFE is clearly able to sharply follow the 
velocity divergence distribution, both on the negative veloc- 
ity divergence side corresponding to infall motions and on 
the positive velocity divergence end representative for the 
outflow from void region. 

It forms a contrast to the rather poor correspondence 
between the density and velocity di vergence distribution i n 
the corresponding reconstruction of IBranchini et all (| 1999T ) : 
the lin-lin plots show that they do not form each oth- 
ers mirror image. We may also note that its velocity di- 
vergence distribution is considerably broader than that of 
the equivalent DTFE pdf. This adheres to a known arte- 
fact of gridbased velocity field interpolations: implicitly 
these yield mass-weighted values whose pdfs have wider 
wings than would be expected for the distribution of proper 
volume-weighted quantities in the mil dly nonlinear regime 
(Bernardeau & van dc WcyEaert 1996). 



Bernardeau fc van de Wevgaertl (|l996l ) convincingly 



demonstrated that the pdf recovered by Delaunay tessel- 
lation (i.e. DTFE) interpolation closely adheres to the pre- 
dictions of Eulerian perturbation theory. The velocity di- 
vergence pdf is strongly sensitive to the underlying cos- 
mological parameters, in particular the cosmological den- 
sity parameter Q m . By applying DTFE towards determin- 
ing the velocity divergence pdf this not only implies the 
possibility to get an accurate estimate of Q m but also 
one which circumvents the usually involved degeneracy be- 
tween the cosmic matter density Q m and the bias b be- 
tween th e matter and galaxy distribution. I n a follow- 
up study (|Bernardeau fc van de WevgaertJll997l ) succesfully 
tested this on a range of N-body simulations of structure 
formation, showing Delaunay interpolation indeed recovered 
the right values for Q m . 

The sensitivity of the pdf to Q m is particularly marked 
near the peak of the pdf and in the maximum void expan- 
sion rate, ie. the maximum value of the velocity divergence 
dmax- The exact location of the peak also includes a depen- 
dence on cosmological scenario and level of nonlinearity. To 
infer f2 m extra assumptions would need to be invoked. It is 
more straightforward to read off the value of Q m from the 
characteristic sharp cut-off of the pdf on the side of the pos- 
itive velocity divergence values. This cut-off relates to the 
maximum expansion rate of voids, which is predicted to be 



= 1.5 fT 



(14) 



in which we ignore the (rather minor) dynamical influence of 
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Figure 14. DTFE velocity shear field projected along the z— supergalactic plane. The slice corresponds to the one presented in Fig. [5] 
We have plotted the velocity shear amplitude a in units of the Hubble parameter. The colour bar indicates the values of the velocity 
shear in the colour contour map. 



a cosmological constant. The value of 1.5 is the difference in 
value of the Hubble parameter in an empty Q — universe 
and that in an Einstein-de Sitter Universe Q — 1, reflecting 
the fact that the interior of the deepest voids locally mimic 
the behaviour of an O = Universe. Following a related idea 
iDekel fc Reel | 19941 ) obtained a strong lower bound on the 
value of Q m ~ 0.3 on the basis of an estimate of the value 
of expansion rate of the nearby Sculptor void. However, for 
a proper determination one should seek to carefully recover 
the sharp cut-off of the pdf. iBernardeau fc van de Weygaertl 



(1996) detailed various arguments why a proper determina- 
tion of the cut-off is nontrivial. They convincingly demon- 
strated the success of tessellation-based methods in outlining 
the pdf's sharp edge. Here we infer from the DTFE deter- 
mined PSCz velocity divergence pdf (fig. [6j that 



fi m = 0.35 ±0.04. 



(15) 



This value agrees with the value determined on the basis 
of the pdf's peak (eqn. 1141 The V ■ v variation within each 
individual bins is used to obtain the error estimate. Given 
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Figure 15. The DTFE velocity shear field on a smoothing scale of 3.8 h — x Mpc (lefthand column) and 10 h _1 Mpc (righthand column). 
The maps in the consecutive rows concern the three mutually perpendicular central supergalactic planes: the X-Y, X-Z and Y-Z plane 
(from top to bottom). Note that the map on a scale of 3.8 h — 1 Mpc is uniform only out to a radius of 30 h~ 1 Mpc, with a gradually 
diminishing resolution towards the outer edge of the sample volume. Also compare the images with the density fields in fig. [8] and fig. [9] 
The greyscale values are in units of the Hubble parameter Hq, with their values given in the bar at the top of the figure. 



that the input velocity field has a value of f3 = 0.5, i.e Q = 
0.315 and 6 = 1, this forms a confirmation of the ability 
of DTFE to self-consistently infer the underlying value of 
cosmological parameters, in particular that of Q m . 



10 VELOCITY SHEAR AND VORTICITY 

The map of the amplitude of the shear component in the 
DTFE velocity field in the supergalactic X-Y plane is shown 
in fig 1141 The map shows the DTFE estimate of the ampli- 
tude of the shear tensor, 

a = E (ffijffii) 1 ' 8 (16) 

where Oij is the symmetric traceless part of the DTFE de- 
termined velocity gradient estimate, 

1 f dvi dvi~] 1 — - N _ ,„„. 

*« = n^+^H^-^- (17) 

Shear flows are induced by the intrinsic asphericity of 
evolving structures and by the external tidal stresses ex- 
erted by the surrounding (inhomogeneous) large scale mat- 
ter distribution. These flows are a manifestation of the 
tidally induced anisotropic contraction of matter into pla- 
nar and filamentary features, locked into a coherent pattern 
through highly dense compact clusters which form at the 
peaks in the primordial density field. It forms the basis fo r 
our understanding of t he Cosmic Web (|Bond et al.lll996l ). 
(van de Wev gaert 200|, also see). 

The impression provided by the shear field in fig. 1141 s 
akin to that of the velocity divergence distribution (fig. Ill[) 
and that of the density distribution (fig[SJ. Evidently there 
is a strong coupling between the density environment and 
the strength of the shear field. The velocity shear attains 
the highest values in and near high density regions, while 
the large empty void regions are regions in which shear 
flows are hardly relevant with respect to the radial out- 
flow (ie. the velocity divergence amplitud e). This is en- 
tirely in line with theoret ical expectations (Hoffman! 1 19861 ; 
iBertschinger fc Jaidfl994h . 

The spatial shear distribution suggests a very strong 
tidal force field in the central region of the Local Universe, 
in the supergalactic X-Y plane. Most outstanding are the 
features in the supergalactic X-Y plane, in particular the 
ridge running through the central Universe from the Pisces 
Perseus supercluster towards the Shapley concentration. On 
the one hand, this impression forms a telling manifestation 
of the characteristic quadrupolar pattern of the local mass 
distribution, defined by the Hydra-Centaurus supercluster 
(the GA region) and the Perseus-Supercluster region. Even 
though undoubtedly real this impression of a strong tidal 
field in the central region is partially biased by the effect 
of a considerably better spatial resolution in the central 
30 /i _1 Mpc of our sample than in the outer regions and the 
resulting absence of the (stronger) higher frequency compo- 



nents present in the center. In order to be able to assess the 
impact of this effect we compare the shear field maps for a 
smoothing radius of 3.8 /i _1 Mpc with those for 12 /i _1 Mpc 
in fig. [TJ] While it confirms the tendency of higher shear 
values in the central region due to the resolved small scale 
contributions, the 12 h~ 1 Mpc convincingly delineates large 
scale patterns responsible for considerable tidal stresses. The 
low resolution, large scale, map does indeed demonstrate 
the reality of the prominent tidal features in the supergalac- 
tic (X-Y)plane. Its prominent dynamical role is additionally 
emphasized by the low resolution X-Z map where the plane 
stands out as the central edge marked by high shear values. 
Interesting features may also be observed in other regions 
of the sample volume: an interesting filamentary extension 
towards [SGY,SGZ]r; [-40, -100] /i _1 Mpc can be observed 
in the Y-Z map. 

10.1 the DTFE Vorticity Map 

Vorticity is not expected to play any role of significance as a 
result of the "linearized" origin of our PSCz velocities. The 
top-hat filter of Rth = 5 /i _1 Mpc used to defined the ve- 
locity sample is an ensurance for the linearity of practically 
all sample velocity data. 

The presence of vorticity in our DTFE reconstructions 
should therefore provide us with a reasonably good impres- 
sion of the influence of systematic artefacts in our maps. 
Regions populated with considerable large or abundant clus- 
ters (e.g. Great Attractor, Pisces-Perseus), do retain some 
measure of vorticity, even after smoothing. Also the DTFE 
method itself also introduces a small spurious vorticity com- 
ponent as a result of its linear interpolation scheme. A care- 
ful analysis of this artefact in our maps, and the effects 
for the DTFE velocity reconstruction, wil l be pr esented in 
iRomano-Diaz. van de Wevgaert fc Schaapl l|2007t) . 



11 SUMMARY AND DISCUSSION 

The Delaunay Tessellation Field Estimator (DTFE, see 
Sch aap fc van de Weygaertl [2000) has been applied to a 
combined analysis of the density and velocity flow field in 
the local Universe. The prime objective of this study has 
been the production of optimal resolution three-dimensional 
maps of the volume-weighted velocity and density fields 
throughout the nearby Universe, the basis for a detailed 
study of the structure and dynamics of the cosmic web at 
each level probed by underlying galaxy sample. In order to 
have a reasonably complete sample of galaxy peculiar ve- 
locities throughout the surrounding local Universe we used 
the PSCz galaxy redshift catalog translated into galaxy po- 
sitions an d velocities by means of the method I linearization 
process of lBranchini et all (|l999T ). 

The spatial galaxy distribution of the PSCz catalog de- 
fines a Voronoi and Delaunay tessellation. These epitomize 



30 Romano-Diaz & van de Weygaert 



the most pure locally defined division of space in the PSCz 
volume. The self-adaptive nature of the Voronoi/Delaunay 
tessellations concerns both spatial resolution and local ge- 
ometry. The high level sensitivi ty to the local point distri- 
butio n is exploited by DTFE dSchaap fc van de Weygaert] 
2000) to produce an estimate of the local density at each 
sample point. The Delaunay tessellation is subsequently 
used as adaptive multidimensional spatial interpolation grid. 
The construction of the Voronoi and Delaunay tessellation 
is the most demanding task of the DTFE routine. Once this 
has been accomplished, they can be used for both the den- 
sity and velocity field (and basically any relevant field). 

Because DTFE is based upon linear interpolation and 
involves constant field gradients within each Delaunay tetra- 
hedron, the analysis of the cosmic velocity field automat- 
ically yields maps of the velocity divergence, shear and 
even vorticity. The spatial distribution of these quantities 
in DTFE is marked by a discontinuous zeroth order map 
marked by Delaunay regions in which divergence, shear or 
vorticity has a constant value. To sensibly relate the DTFE 
velocity divergence and shear fields it is necessary to smooth 
over a scale comparable to that of the corresponding Delau- 
nay cells. 

Earlier studies l|Schaap fc van de Weygaert] [2006 ) have 
made clear that DTFE is particularly optimalized for the 
analysis of mass distributions marked by one or more of the 
essential aspects of the nonlinear cosmic web of the Mega- 
parsec Universe (|Bond et al.| [l996). One of the main advan- 
tages of DTFE is its large dynamic range, allowing it to 
trace small scale structures along with the large scale en- 
vironment in which they are embedded. The spatial and 
morphological adaptivity of DTFE is ideally suited for dis- 
secting the essential components of the nonlinear weblike 
cosmic matter distribution. It resolves a hierarchically struc- 
tured matter distribution to the smallest possible resolution 
scale set by the particle number density. Perhaps most out- 
standing is the ability of DTFE to retain the morphology 
cq. shape of the features and patterns in the matter distri- 
bution: the characteristic anisotropic filamentary and pla- 
nar features of the Cosmic web are fully reproduced in the 
conti nuous DTFE density field (|Schaap fc van de Weygaert] 
2000). Also the third major characteristic of the nonlinear 
Megaparsec universe, the dominant existence of near-empty 
voids, is resolved by DTFE. Both their flat internal density 
distribution as well as their sharp outline and boundary are 
recovered in detail, while shot noise in these sparsely sam- 
pled regions tends to be suppressed. 

The maps of the PSCz local Universe reveal a sharply 
defined density field in which many familiar structures and 
features can be recognized at an optimal spatial resolu- 
tion. The region around the Local Group and its filamen- 
tary extension towards the Hydra-Centaurus and Shapley 
complexes as well as the Pisces Perseus supercluster are 
prominently resolved in the DTFE density fields. Moreover, 
these features are as clearly and sharply defined in the corre- 
sponding velocity flows. Perhaps the most outstanding fea- 
tures in the velocity vector fields are the radially expand- 
ing void regions. In comparison to earlier published maps of 
the d ensity field and veloc i ty flow in the PSCz sample vol- 
ume (Bran chini et al ] |l999l ; ISchmoidt" et al . 1999) the DTFE 
maps have a considerably sharper defined appearance. While 
conventional maps suppressed small scale details and were 



at loss in undersampled void regions, in particular wrt the 
interpolation of the velocity field, the DTFE fields reveal a 
beautifully textured flow field. The DTFE flows are marked 
by prominent bulk flows, shear flows and radial inflow (Coma 
cluster) . Arguably the most outstanding and unique feature 
of the DTFE maps is the sharply defined radial outflow re- 
gions in and around underdense voids, marking the dynam- 
ical importance of voids in the Local Universe. 

The present study includes two specific aspects and 
challenges for DTFE: (1) the application of DTFE to recon- 
struct the density field as well as the corresponding velocity 
field in the same cosmic volume and (2) the gradually dimin- 
ishing spatial resolution as a function of radial distance as a 
consequence of the flux-limited nature of the PSCz sample. 

While one may correct estimated density values by tak- 
ing into account the well-defined radial selection function, it 
is not possible to correct for the loss in spatial resolution: 
the lower galaxy sampling density goes along with a loss, 
at increasing distances, in spectral coverage of the maps. 
A full recovery would have to involve pre-conceived notions 
on the cosmological structure formation scenario, cq. the 
power spectrum. The relatively prominent central values in 
the density, velocity divergence and shear maps are a re- 
flection of this effect. On the other hand, the velocity field 
itself - less affected by the lack of high frequency spectral 
power at large distances - for all practical purposes appears 
to be uniformly covered by the DTFE reconstruction. The 
sharply outlined edges of superclusters and the radial out- 
flow of voids are found throughout the PSCz volume. 

We have addressed the relationship between the den- 
sity field and velocity field by correlating the density values 
with the local velocity divergence and shear. Overall, the 
linear density-velocity divergence relationship is accurately 
reproduced. This reflects the origin of the galaxy velocities, 
having been computed from a redshift map through a lin- 
earization procedure. These simple quantitative relations do 
not express the spatial coherence and correlations within 
the velocity field. The tight spatial correspondence between 
density features and velocity flows in the DTFE maps - 
voids with radial outflow, elongated and flattened super- 
clusters with bulk and shear flows - does demonstrate the 
remarkable adaptive nature of the DTFE technique and its 
promise of understanding the dynamics of the cosmic web. 
The maximum expansion rate of voids defines a sharp cutoff 
in the velocity divergence pdf. It enabled us to test the self- 
consistency of the DTFE method. Indeed it confirmed the 
value of £l m ~ 0.35 of the cosmology underlying the galaxy 
velocity data. 

For furthering our insight into the velocity field in the 
nearby universe it would be desirable to be able to probe 
quasi-linear contri butions to the velocity field. Techniques 
like the FAM— z of iBranchini et al.l (|2002f) . applied in com- 
bination with the DTFE algorithm, would offer a potentially 
promising approach. While this study has shown the poten- 
tial of DTFE towards a succesfull analysis of velocity fields 
it concerns a test on the basis of a sample of pre-processed 
galaxy velocities. This in cludes the assumption that the ve- 
locity field linearization (IBranchini et al.l 1 1999T ) has succes- 
fully dealt with redshift-space distortions and sampling ef- 
fects. 

Ultimately we will seek to develop a DTFE-based for- 
malism to directly analyze samples of galaxy peculiar ve- 
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locities on the basis of the "raw" observational data. The 
challenge is found in a few aspects. One aspect concerns er- 
rors and sampling effects. A second aspect involves redshift 
distortions. A third complicating aspect involves fundamen- 
tal assumptions behind the velocity field reconstruction pro- 
cedure, such as the requirement of laminarity of the flow. 
An extensive study of various sampling effects, such as de- 
creased sampling density and inhomogeneous sampling, on 
simulation-based data do show that the noise level of the 
DTFE velocity reconstructions does i ncrease but that thi s 
does not involve any systematic shifts (I Romano- Dfadl2004 ). 
Given the acclaimed optimal triangulation properties of De- 
launay grids underlying the development of DTFE this may 
not come as a surprise. The same study involved tests of 
the effects of a variety of different error contributions. They 
induced similar effects as a less balanced galaxy sampling. 
These results will be presented in an extensive study and 
report on error and sampling effects in a forthcoming paper. 
More challenging issues concern the physically more pro- 
found aspects of redshift distortions and higher-order flow 
characteristics. As a first-order version o f the wider c l ass of 
natural- neig hbour interpolation schemes ((Sibson f 980. fl98ll ; 
IWatsorj|l992l ). DTFE may be extended to a higher-order 
natural neighbour interpolation scheme based upon the De- 
launay tessellation of a galaxy sample. Such a formalism 
would enable a more complex treatment of nonlinear velocity 
flows, possibly even of multistream flows. Already we have 
been working towards the impl ementation of these higher- 
order interpolation schemes (see lvan de Weygaert fc Schaad 
120071 . for a review on the basics), and we plan to look into the 
possibility on the basis of routines made available through 
CGAL 7 . On the other hand, highly mixed virialized regions 
will remain beyond the scope of these interpolation tech- 
niques. This may hardly be considered a constraint given 
the fact that for the foreseeable future the focus of nearly all 
galaxy peculiar velocity surveys will be the dynamics of large 
scale structures. The inversion from redshift to real space 
data may be facilitated and enabled by using the virtues of 
a Delaunay grid for the solution of PDEs. The use of Delau- 
nay tessellations as grids f or the numerical solution of PDEs 
has been first described bv lBraun fc Sambridgel(|l995h . their 
success provides substantial imperative for the formulation 
of the redshift space to real space inversion following these 
lines. 

In the meantime DTFE has been elaborated upon 
through the development of a num ber of specific fea - 
ture detecting techniques. The MMF l|Aragon et all 120071 ) 
forms a highly promising tool for selection and identi- 
fication of filamentary and planar structures and their 
spatial relationship with neighbouring structures. The 
watershed based void detection algori thm developed by 
iPlaten. van de Weygaert fc Jonesl (2007) has been applied 
tow ards the analysis of the hierarch ical evolution of voids 
(see ISheth fc van de Weygaert! |2004| ). In all cases, the fully 
adaptive and optimal tessellation based characteristics of 
DTFE form the crucial starting point. 



7 CGAL is a C++ library of algorithms and data structures for Com- 
putational Geometry, see www . cgal . org 
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